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We present and compare projected distributions of mass, galaxies, and the intrachister medium (ICM) 
for a sample of merging clusters of galaxies based on the joint weak-lensing, optical photometric, and X- 
ray analysis. Our sample comprises seven nearby Abell clusters, for which we have conducted systematic, 
' deep imaging observations with Suprime-Cam on Subaru telescope. Our seven target clusters, representing 

various merging stages and conditions, allow us to investigate in details the physical interplay between dark 
matter, ICM, and galaxies associated with cluster formation and evolution. A1750 and A1758 are binary 
systems consisting of two cluster-sized components, A520, A754, A1758N, A1758S, and A1914 are on-going 
Q..^ ' cluster mergers, and A2034 and A2142 are cold-front clusters. In the binary clusters, the projected mass, 

I optical light, and X-ray distributions are overall similar and regular without significant substructures. 

\^ I On-going and cold-front merging clusters, on the other hand, reveal highly irregular mass distributions. 

' Overall the mass distribution appears to be similar to the galaxy luminosity distribution, whereas their 

, distributions are quite different from the ICM distribution in a various ways. We also measured for 

individual targets the global cluster parameters such as the cluster mass, the mass-to-light ratio, and the 
ICM temperature. A comparison of the ICM and virial temperatures of merging clusters from X-ray and 

■ weak-lensing analyses, respectively, shows that the ICM temperature of on-going and cold-front clusters is 
significantly higher than the cluster virial temperature by a factor of ~ 2. This temperature excess in the 

Q ■ ICM could be explained by the effects of merger boosts. 

5-H ' Key words: gravitational lensing — X-rays: galaxies: clusters — galaxies: clusters: individual (A520, 

to '. A754, A1750, A1758, A1914, A2034, A2142) 

- ^ , -. INTRODUCTION ter. The merger-induced shocks and turbulences generate 

■ relativistic particles, which then produce synchrotron ra- 
I In hierarchical structure formation models based on the dio emission and inverse-Compton emission on the Cosmic 

Cold Dark Matter (CDM) scenario, small structures form Microwave Background. We however have not yet under- 
first, and they merge to form larger, more massive ob- stood the cluster merger phenomena and the relationships 
jects. Since clusters of galaxies are the largest gravita- between the ICM and dark matter in detail, 
tionally bound objects in the universe, they are still form- Observationally, cluster mergers have been investi- 
ing through mergers of sub-clusters and smaller groups gated based on optical measures of cluster dynamics, 
of galaxies along the filamentary structure within which namely line-of-sight velocity distributions of cluster mem- 
they are embedded. The cluster merger is one of the ber galaxies as well as optical appearance of substruc- 
most energetic events in the universe because of its high- tures. Optical evidence of substructures were reported by 
est binding energy and huge energy release. Dark mat- many authors based on apparent concentrations of galax- 
ter, which is the dominant mass component of clusters, ies in projected space (e.g., Abell, Neyman, & Scott 1964; 
governs the dynamics of cluster mergers, and drives var- Baier & Ziener 1977; Geller & Beers 1982), although a 
ious phenomena in the hot, diffuse intracluster medium detection of subclustering of cluster galaxies is not direct 
(ICM). Colliding substructures with high collision veloc- evidence for the existence of underlying mass component, 
ities trigger shocks or turbulences in the ICM. Merger On the other hand, one can investigate the galaxy dynam- 
shocks increase the temperature and entropy of the ICM, ics in clusters by examining distribution of galaxies both 
resulting in increasing the X-ray luminosity of the clus- in spatial and velocity space (e.g. Zabludoff et al. 1993; 

Ashman, Bird & Zepf 1994), 



Based on data collected at Subaru Telescope, which is operated 
by the National Astronomical Observatory of Japan. 



In recent years. X-ray satellites such as ASCA, ROSAT, 
Chandra, XMM-Newton and Suzaku make it possible to 
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study the physical processes in the ICM involved in the 
formation and evolution of clusters. In particular, precise 
spectro-imaging data of Chandra and XMM-Newton al- 
low us to derive the temperature, pressure, and entropy 
maps of the ICM, which are useful tools to study the clus- 
ter merger process. Recent Chandra and XMM-Newton 
observations have revealed a highly-complex X-ray mor- 
phology of merging clusters (e.g. Govoni et al. 2004; 
Henry, Finoguenov, & Briel 2004; Finoguenov, Bohringer, 
& Zhang 2005). Markevitch et al. (2000) discovered 
an unexpected ICM structure of contact discontinuities, 
which are called "cold fronts". However, although such 
X-ray observations of the ICM physics provide us with 
detailed information of the cluster merger, X-ray study 
alone is not sufficient for understanding of complex phys- 
ical processes acting during the cluster merger. This is 
because cluster mergers are dynamically governed by dark 
matter which mainly accounts for ^ 80% of the total mass 
of clusters. Further, the ICM in merging clusters is not 
in, or close to, hydrostatic equilibrium, so that the X-ray 
determined total mass of merging clusters, under the as- 
sumption of hydrostatic equilibrium, will be different from 
the gravitational mass of the clusters. 

Gravitational lensing effects on background galaxies, 
on the other hand, are powerful, unique tools to study 
the mass distributions in clusters regardless of the physi- 
cal/dynamical state of matter in the systems. Weak lens- 
ing image distortions of background galaxy images can be 
used to map the distribution of matter in clusters on scales 
of ~ 100 kpc to 1-2 Mpc scales in a model-independent 
way (see, e.g., Bartelmann & Schneider 2001; Umetsu, 
Tada, & Futamase 1999). Therefore weak lensing enables 
the direct study of mass in clusters even when the clusters 
are in the process of pre/mid/post merging, where the as- 
sumption of hydrostatic equilibrium or isothermality are 
no more valid. 

In this way, weak lensing. X-ray, and optical photomet- 
ric observations provide complementary information of the 
physical properties of clusters. Hence, a joint analysis of 
weak lensing. X-ray, and optical photometric observations 
will yield a comprehensive picture of the physical state of 
ICM. dark matter and member galaxies. 

In this paper we present weak-lensing mass distribu- 
tions of a sample of seven nearby merging clusters, re- 
constructed using weak shear data taken with Suprime- 
Cam on the 8.2m Subaru telescope, and make a quan- 
titative comparison of the mass maps with X-ray and 
optical galaxy distributions. The Subaru/Suprime-Cam 
is the most ideal, working instrument for weak lensing 
shape measurements of background galaxies because of its 
wide field-of-view of 34' x 27', photon collecting power of 
the 8.2m mirror, and excellent image quality with stable, 
small PSF anisotropy (e.g., Miyazaki et al. 2002; Hamana 
et al. 2003; Sato et al. 2003; Broadhurst et al. 2005). 
We selected for the Subaru weak lensing observations 
seven nearby merging clusters (0.05 ^ z ^ 0.28) of differ- 
ent merging stages and properties. A520, A754, A1750; 
A1758, A1914, A2034 and A2142, for which archival X- 
ray data of Chandra and/or XMM-Newton are available. 



Detailed analysis results of individual clusters will be pre- 
sented elsewhere. Details of our X-ray, optical and weak- 
lensing analyses are presented in §2. In §3, we compare the 
resulting projected distributions of mass, optical luminos- 
ity density, and X-ray emission for our sample of merging 
clusters. In §4 we derive global cluster parameters such as 
the cluster mass, the cluster mass-to-light ratio, and the 
ICM temperature. Finally, in §5, we summarize the pro- 
jected offset distributions of cluster galaxies and ICM with 
respect to the mass, and discusss the implications of our 
results for understanding of cluster merger phenomena. 

Throughout the paper we adopt the concordance cos- 
mology of n,no = 0.3, flA = 0.7 and = 70 kms"^Mpc"\ 
and use the AB magnitude system. 

2. DATA ANALYSIS 

2.1. Data Preparation 

2.1.1. Subaru Data Analysis 

We observed all of the target clusters, except A520 
for which archival data are available (SMOKA), with 
Subaru/Suprime-Cam (ID: S05A-159, PI: N. Okabe). For 
each target cluster we chose Rc or i' band filters for the 
weak lensing shape analysis, and g' or V for separat- 
ing member and background galaxies based on the color- 
magnitude diagram. All of our observations were taken 
using the AG (acquisition and guide) probe for guide stars, 
while the archival i'-band data of A520 were taken with 
and without guide probe on 17th November 2001 and 19th 
October 2001, respectively. We retrieved from SMOKA 
a total of seven i'-band images taken under good seeing 
conditions (~ 0."6), and analyzed separately the following 
three imaging data sets: (A) all of seven i' images taken 
both with and without guide probe (7 x 240s exposure), 
(B) four i' images taken without guide probe (4 x 240s ex- 
posure), and (C) three i' images taken with guide probe 
(3 X 240s exposure). Our main analysis results are based 
on the deepest data set (A) (see §3.7), which was also 
used by Mahdavi et al. (2007) for their multi-telescope, 
multi-bandpass weak lensing analysis of the cluster. We 
compare in Appendix the resulting mass maps derived 
from the thee different imaging data sets. We note that 
the A520 data were taken with large dithering offsets of 
« 2.'3, whereas the other cluster data were taken with a 
dithering offset of 1'. A summary of observation parame- 
ters and conditions is given in Table 1. 

Each Subaru data set was reduced using the software 
SDFRED (Yagi et al. 2002; Ouchi et al. 2004). For each 
CCD frame we first estimated and subtracted the bias 
level by the median value in the overscan region. Next 
we carried out flat-fielding for all of the CCD frames to 
correct for the pixel-to-pixel variation in sensitivity. We 
then combined mosaic CCD images and corrected for the 
field distortion due to the camera optics. In order to ob- 
tain the local sky background map we first measure the 
median sky level in each mesh of a grid, and then compute 
a high-resolution background map via the 2D bilinear in- 
terpolation. Here we took the mesh size to be 32 x 32 pix- 
els with 0"202 pixel"^. We subtract the sky background 
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measured from our images. We matched the point spread 
function (PSF) between the frames and masked out areas 
and objects vignetted by the AG probe and bad pixels, 
such as satellite trails. Finally we generate a median com- 
bined image from these individual frames. Here stacking 
parameters such as the coordinate shifts, rotations, and 
scalings were determined by stellar objects common to all 
exposures. A summary of the seeing FWHMs of final re- 
duced images for the weak lensing analysis is shown in 
Table 1. 

In order to make an accurate comparison between op- 
tical and X-ray images, we performed astrometry using 
the IRAF 2.11 tasks (CCMAP and CCSETWCS) with 
our object catalogs generated using SExtractor (Berlin & 
Arnouts 1996) and the USNO-A2.0 catalog as our refer- 
ence. After the astrometry correction typical rms residu- 
als are about 1 pixel (=0."202). 

2.1.2. Chandra Data Analysis 

Wc used archival Chandra ACTS data for the X-ray 
analyses of A520, A754, A1914 and A2142. The X-ray ob- 
servations are identified by their ID numbers (Obs. ID) in 
Table 1. The data were reduced using the CIAO (version 
3.3 for A520, A754 and A2142, and ver 3.2.2 for A1914) 
for the X-ray data processing. Standard screening was 
applied to the photon list. Bad pixels and columns were 
removed with ASCA grades 0, 2, 3, 4 and 6. The data 
were then cleaned of periods of anomalous background 
rates which are more than 3ct significance from the mean 
value. We found an aspect offset in the reduced data of 
A754 (Obs. ID = 577), which was then corrected using the 
aspect calculator. We applied the CTI correction to our 
data with focal plane temperature of — 120°C, because of 
the charge loss in ACIS-I CCD chips during transfer to 
the readout node. Background was derived from blank 
sky data. Wc normalized the background data so that the 
count rates in the 10— 12 keV band matched our observed 
count rates in the outer region of Chandra ficld-of- views. 

The raw X-ray image was binned by a factor of 8 in each 
dimension (DETX and DETY). Point sources detected by 
WAVDETECT in the CIAO package was excluded from 
the images. X-ray images were adaptively smoothed, and 
background and exposure images were smoothed using the 
same kernel. Smoothed X-ray images were obtained by 
substituting the smoothed backgrounds and dividing them 
by smoothed exposure maps. 

2.1.3. XMM-Newton Data Analysis 

Wc used archival XMM-Newton data for A1750, A1758 
and A2034 from three EPIC cameras (MOSl, M0S2 and 
PN). X-ray observations are identified by their ID num- 
bers (Obs. ID) in Table 1. We created calibrated event 
files using SAS version 6.5.0. The XMM-Newton data 
during high background flares were removed. To identify 
the high background periods, we generated light curves 
in the > 10 keV energy band. The light curves were 
binned in 100 s intervals. We obtained good time inter- 
vals (GTIs) by applying 2a clipping. Wc confirmed that 
GTIs are consistent with previous works done by Bclsolc 
etal. (2004) and David &Kempncr (2004). The data were 
then filtered to leave only events with PATTERN < 12 and 
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#XMMEA_EM for MOS and FLAG = and PATTERN < 
4 for PN. 

We followed the double-subtraction method (Arnaud 
et al. 2002) for background modeling. Background files 
were produced by taking into account three different back- 
ground components: particle-induced background, cosmic 
X-ray background (CXB) and residual soft photon con- 
tamination. The particle background was estimated from 
the filter wheel closed (FWC) data released from the EPIC 
Background Working Group. The particle background 
was re-normalized in the 0.3 — 10 keV energy band by 
the ratio of the observation counts over the FWC particle 
counts outside the XMM field-of-view. For a CXB back- 
ground estimation, we used blank-sky background data 
of the same instrument, filter and mode as those used 
in each observation. The blank-sky background counts 
were normalized to match the observed event counts in 
a source-free annulus at large radius. We obtained soft 
X-ray excess images by matching in the low energy band 
the background count rates to our observed count rates in 
the outer region of the XMM field-of-view. Finally these 
normalized background images were combined to form the 
total X-ray background. 

The raw X-ray image was then smoothed by a factor of 
64 in each dimension (DETX and DETY) . Finally wc ob- 
tained for each target adaptively smoothed, background- 
subtracted and exposure-corrected mosaic of the MOSl, 
M0S2 and PN images, in the same manner as done for 
the Chandra data. 

2.2. Cluster Galaxy Selection 

We performed object detection and aperture photome- 
try using SExtractor (Berlin & Arnauts 1996) in dual- 
image mode with deep i' or Rc as the detection im- 
age, yielding a common object catalog for the two filters, 
{g',Rc) or {V,i'). We extracted all objects with isopho- 
tal area larger than 10 pixels above 3a pixel"^ of the lo- 
cal sky level. We used MAG.AUTO and MAG.APER of 
SExtractor output as object's total magnitude and aper- 
ture magnitude, respectively. For aperture photometry 
the aperture size is set to 10 pixels (2."02), except 20 pix- 
els (4"04) for A754 at a very low redshift of z = 0.0542. 
We define for each target cluster a sample of cluster mem- 
ber galaxies according to their color {g' — Rc or V — i') and 
magnitude {Rc or i'). Figure 1 shows the color magnitude 
diagram (CMD) for one of our target clusters, A2142, af- 
ter the removal of stellar objects. The CMD in Figure 1 
exhibits a tight color-magnitude (CM) sequence of early- 
type cluster galaxies, since the two filters bracket the rest- 
frame 4000A break. We quantify the CM sequence of 
galaxies with a linear relation in color-magnitude space. 
We set the maximum magnitude of the linear CM rela- 
tion to Rc,i' = 22 ABmag. It can be seen that the lumi- 
nous end of CM sequence extends to the high-luminosity 
regime of Rc « 16 ABmag. For individual galaxies, we 
convert from apparent to absolute magnitudes by using 
the k-corrcction for early-type galaxies. Here we assume 
that all of the member galaxies in each target cluster are 
located at a single redshift summarized in Table 2. 
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We measure the optical luminosity of an individual clus- 
ter target by summing all the cluster galaxies within a 
given radius from the cluster center. Here we define the 
cluster center to be either at the brightest cluster galaxy 
(hereafter BCG). To do this we use a sample of bright CM- 
sequence galaxies with magnitudes Rc,i' < 22 ABmag. We 
note that BCGs tend to deviate from the linear CM re- 
lation. Hence we visually checked if such brightest galax- 
ies were properly included in our cluster member sample, 
and included them if they were missing. A field correc- 
tion is applied to the measured luminosity density of clus- 
ter galaxies to account for contamination by field galax- 
ies: We estimate the background luminosity density from 
an annular region outside of the target region, and then 
subtract this background contribution from the observed 
luminosity density of the cluster. 

In order to measure the total luminosity of cluster 
galaxies, we also need to correct for incompleteness of 
the sample due to the magnitude limit. We assume that 
cluster galaxies follow a Schechter luminosity function 
(Schechter 1976) of the form: 



(1) 



By integrating (l){L) down to the luminosity cutoff corre- 
sponding to the magnitude limit Lmin of our cluster mem- 
ber sample, we obtain the following relation between the 
total luminosity Ltot and the observed luminosity Lobs of 
the cluster: 



it 



r(2-p) 



T(2-p,Lum/i*)^ 



(2) 



where T{a,x) is the upper incomplete gamma function. 
We adopt the following parameters for the cluster lumi- 
nosity function: p^im,M^ = -21.89 + 51og(/i) in R^- 
band and p = 0.70, M*, = -22.31 + 5 log(/i7o) in i'-band 
(Goto et al. 2002), where M* is the absolute magni- 
tude corresponding to the characteristic luminosity L*. 
The resulting luminosity correction factors, T{2~p) /r(2 — 
P, Liim/ L*), are of order of unity as shown in Table 3. 

2.3. Weak Lensing 

In the weak lensing analysis we aim to reconstruct from 
image distortions of background galaxies the dimension- 
less surface mass density 

A^(0) = S-/I](0) (3) 

of clusters in units of the critical surface mass density for 
gravitational lensing, 

^rr ^ . „ ^ 1 . (4) 



47rG DdDas ' 

where Ds , Dd and Dds are the angular diameter distances 
from the observer to the sources, from the observer to 
the deflecting lens, and from the lens to the sources. In 
the following we closely follow the standard notation of 
Bartelmann & Schneider (2001). 



2.3.1. Weak Lensing Distortion Analysis 

We use for our weak lensing analysis the IMCAT pack- 
age developed by N. Kaiser ^, following the formalism de- 
veloped by Kaiser, Squires, & Broadhurst (1995). We 
have modified the method somewhat following the proce- 
dures described in Erben et al. (2001). We used the same 
analysis pipeline as in Broadhurst, Takada, Umetsu et al. 
(2005), Medezinski, Broadhurst, Umetsu et al. (2007), 
and Umetsu & Broadhurst (2007). Note that our weak 
lensing analysis on A520 is based on the deepest data set 
in Table 2, and a comparison of the mass reconstructions 
between different imaging data sets is given in Appendix, 
demonstrating the effects of different observation condi- 
tions and strategies on the weak lensing analysis. 

We measure the image ellipticity = 
{Qii - Q22,Qi2} /(Qii + Q22) from the weighted 
quadrupole moments of the surface brightness of individ- 
ual galaxies, 
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d^ew{e)eaOpi{e) (a,/3 = i,2) 



(5) 



where I{9) is the surface brightness distribution of an ob- 
ject, W{9) is a Gaussian window function matched to the 
size of the object. 

Firstly the PSF anisotropy needs to be corrected using 
the star images as references: 



pQ/3 * 



(6) 



where Psm is the smear polarisability tensor being close 
to diagonal, and g* = {P*sm)af!^* is the stellar anisotropy 
kernel. We select bright, unsaturated foreground stars 
identified in a branch of the circularized half-light radius 
(r/i) vs. magnitude {i') diagram to measure . In Table 2 
we summarize basic information and statistics of the stel- 
lar samples used for PSF corrections in our weak lensing 
analysis. To obtain a smooth map of g* used in equation 
(6), we divided the co-added mosaic image, whose size is 
about ^ llK X 9K pixels, into 5x4 chunks. The chunk 
length is determined based on the typical coherent scale of 
PSF anisotropy patterns. By this, PSF anisotropy in in- 
dividual chunks can be well described by fairly low-order 
polynomials. The typical number of stars per chunk is 
20-50, depending on the field (see N* of Table 2). We 
fitted the q* in each chunk independently with second- 
order bi-polynomials. q'^{0), in conjunction with iterative 
cr-clipping rejection on each component of the residual: 
5e* = e* — P*fmq*p{(i)- As summarized in Table 4, un- 
corrected ellipticity components of stellar objects have on 
average a mean of 1 — 2% with a few % of rms, or variation 
of PSF across the data field. On the other hand, the resid- 
ual (5e* after correction is reduced to |<5e^| ^ 10~^. After 
the anisotropic PSF c orrectio n the rms value of stellar el- 
lipticities, u(5e*) = \/ (|(5e*|^), is reduced from a few % to 
(4 — 8) X 10~^. We show in Figure 3 distributions of stellar 
ellipticity components before and after the PSF anisotropy 
correction for our target clusters. Figure 4 shows the dis- 
tortion fields of stellar ellipticities before and after the 
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PSF anisotropy correction. From the rest of the object 
catalog, we select objects with nT* ^ rh < 15 pixels as 
a weak lensing galaxy sample, where nT* is the median 
value of stellar half-light radii (see Table 4), correspond- 
ing to the half median width of circularized PSF over the 
data field. An apparent magnitude cut off is also made to 
remove from the weak lensing galaxy sample bright fore- 
ground/cluster galaxies and very faint galaxies with noisy 
shape measurements. Table 5 summarizes the magnitude 
range and the mean surface number density Ug of back- 
ground galaxies for our sample of target clusters. Without 
color selection the Ug is ranging from 7ig « 37arcmin^^ 
(A520) to Ug w 72arcmin"^ (A2142), depending on the 
depth, the seeing condition of the observations, and the 
degree of contamination of cluster/foreground galaxies. 

Second, we need to correct the isotropic smearing effect 
on image ellipticities caused by seeing and the Gaussian 
window function used for the shape measurements. The 
reduced shear = 7^/(1 — k) can be estimated from 

9a = {Pg~')a/3e'p (7) 

with the pre-seeing shear polarizability tensor P^p- In the 
weak lensing limit where |At|, |7q| <d, the image ellipticity 
is linearly proportional to the gravitational shear 70,: ga~ 
la ~ {Pg^)a0e-'0- We follow the procedure described in 
Erben et al. (2001) to measure Pg for an individual galaxy 
(see also § 3.4 of Hetterscheidt et al. 2007). We adopt the 
scalar correction scheme, namely, 
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-tr[Pg]Sc,p = PIS^^ 



(8) 



(Hudson et al. 1998; Hoekstra et al. 1998; Erben et 
al. 2001; Hettersheidt et al. 2007). However, the 
measured for individual galaxies are still noisy especially 
for small and faint galaxies. We thus adopt a smooth- 
ing scheme in object parameter space proposed by Van 
Waerbeke et al. (2000; see also Erben et al. 2001; Hamana 
et al. 2003). In this scheme we first identify N neighbors 
for each object in r^-mag parameter space. We then cal- 
culate over the local ensemble the median value {PI) of 
Pg and the variance ti^ of 
The dispersion dg 
estimate for individual galaxies. We adopt = 30. For 
each cluster field we compute the mean variance ct^ over 
the background galaxy sample. In Table 5 we listed the 

mean rms Og = per galaxy, which is of the order of 
ffg ~ 0.4. Finally we use the following estimator for the 



g 9 ~ 9i+ 192 using equation (7) . 
is used as an rms error of the shear 



reduced shear: 



'al{Pl)^ 



(9) 



9a = 

2.3.2. Weak Lensing Mass Reconstruction 

Having obtained the shear estimates for a sample of 
background galaxies, we then pixelize the distortion data 
into a regular grid of pixels using a Gaussian 'Wg{9) oc 
exp[-6'V6'g] with Og = FWHM/V41n2. Further we incor- 
porate in the pixelization a statistical weight Ug for an 
individual galaxy, so that the smoothed estimate of the 
reduced shear field at an angular position 6 is written as 



Y.^^g((^-e^) 



(10) 



where ga.i is the reduced shear estimate of the ith galaxy 
at angular position 0^, and Ug_i is the statistical weight 
of ith galaxy which is taken as the inverse variance, 
Ug,i ~ ^/i'^g i +0!^); with (Tg^i bclug the rms error for the 
shear estimate of ith galaxy (see § 2.3.1) and being the 
softening constant variance (Hamana et al. 2003). We 
choose a = 0.4, which is a typical value of the mean rms 
(Tg over the background sample. The case with a = 
corresponds to an inverse-variance weighting. On the 
other hand, the limit a 3> CFg^i yields a uniform weight- 
ing. We have confirmed that our results are insensitive to 
the choice of a (i.e., inverse- variance or uniform weight- 
ing) under the adopted smoothing parameters. 

The error variance for the smoothed shear g — gi + ig2 
(10) is then given as 



(.T,^Wg.^Ug,^) 



(11) 



where Wg,i = Wgid — 6i) and we have used {ga,ig/3,j) = 
{l/2)aU^/^Sfj with S^^^ and Sf^ being the Kronecker's 
delta. 

We then invert the pixelized reduced-shear field (10) to 
obtain the lensing convergence field. In the map-making 
we assume the linear shearing in the weak- lensing limit, 
that is, ga = Ja/ (1 — k) « 7^. We adopt the two inversion 
methods, namely, the Kaiser & Squires inversion method 
(Kaiser & Squires 1993) and the noise-filtering inversion 
method outlined in Seitz & Schneider (2001). The first 
method makes use of the 2D Green function in an in- 
finite space, while the latter is based on the finite-field 
solution of the inversion problem for the reconstruction 
kernel. The finite-field method must be used for a nearby 
cluster where the data field is dominated by positive, bi- 
ased density field with (k) > 0. In the linear map-making 
process, the pixelized shear field is weighted by the inverse 
of the variance (11). Note that this weighting scheme 
corresponds to using only the diagonal part of the noise 
covariance matrix, N{9i,dj) ~ {Ag{di)Ag{Oj)), which is 
only an approximation of the actual inverse noise weight- 
ing in the presence of pixel-to-pixel correlation due to non- 
local Gaussian smoothing. In Table 5 we summarize the 
Gaussian FWHM used in the pixelization and the rms 
noise level in the reconstructed k field for our sample of 
target clusters. 

The smoothing scale is chosen so as to optimize the 
weak lensing detection of target mass structures, depend- 
ing both on the size of the structure and the strength 
of noise power (cx ag/rig). Under the adopted smoothing 
parameters, the typical rms reconstruction error in the k 
map is ~ 0.02 except for A1914, in which a subarcmin- 
scale (FWHM = 0f75), high resolution reconstruction was 
adopted to resolve substructures (see §3.4). 
2.3.3. Red Background Galaxy sample 

As demonstrated by Broadhurst et al. (2005) and 
Medezinski et al. (2007), it is crucial to make a secure 
selection of background galaxies in order to minimize the 
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dilution of the Icnsing signal by cluster/foreground galax- 
ies and to make an accurate determination of the cluster 
mass. To do this, wc define a sample of red background 
galaxies whose colors arc redder than the CM sequence of 
cluster member galaxies due to large /c-corrections. These 
red background galaxies are largely composed of early to 
mid-type galaxies at moderate redshifts (Medezinski et al. 
2007). Cluster member galaxies are not expected to ex- 
tend to these colors in any significant numbers because 
the intrinsically reddest class of cluster galaxies, i.e. E/SO 
galaxies, are defined by the CM sequence and lie blueward 
of chosen sample limit, so that even large photometric er- 
rors will not carry them into our red sample. 

In Table 6 we summarize our selection criteria for the 
red background sample and the resulting mean surface 
density of the red galaxies. The larger the cluster red- 
shift, the redder the CM sequence of cluster galaxies due 
to larger fc-corrections. Hence, the number of red back- 
ground galaxies, selected in this way, will decrease with 
cluster redshift. In particular, A1758 is at a moderately 
high redshift of z = 0.279, and there are very few such red 
galaxies remained. We therefore relaxed the selection cri- 
teria of the red background sample for the case of A1758 
(see Table 6). 

2.3.4- Tangential Distortion and Cluster Mass Profile 

For an individual cluster, we derive an azimuthally- 
averaged shear profile as a function of projected radius 
from the fiducial cluster center, which is chosen as the 
optical center (see §2.2). 

The tangential component is used to obtain the az- 
imuthally averaged distortion due to lensing, and com- 
puted from the distortion coefficients (51,(72) of each ob- 
ject: 

.9+ = -(gicos20-|-g2sin2(/)), (12) 

where is the position angle of an object with respect to 
the cluster center, and the uncertainty in the 5+ measure- 
ment is (T = <7g/\/2 in terms of the rms error ag for the 
reduced shear measurement. 

The estimation of 5+ only has significance when evalu- 
ated statistically over a large number of background galax- 
ies because of the intrinsic spread in shapes and orienta- 
tions as well as the measurement errors in the shape mea- 
surement. In radial bins we calculate the weighted average 
of the g^s and its weighted error: 

(.9+(g„)) = (13) 

CT+(en) = l/^H^^, (14) 

where the index i runs over all of the objects located 
within the nth annulus with a median radius of 9m and 
""g.i = i + '^^) inverse variance weight softened 

with a = 0.4 (sec §2.3.2). 

For a parameter-free estimation of the cluster mass pro- 
file we use the aperture-densitometry, or so-called the Cc- 
statistic (Fahlman et al. 1994; Clowe et al. 2000) of the 
form: 



Cc(e;^i„n,^out) - 't(< e) - RiOinn < < ^out) (15) 

= 2 [ '"°dln(?'(7+(0)) 
Jo 

where 0inn and ^out are the inner and outer radii of the 
annular background region in which the mean background 
contribution, K(0inn <0 < ^out)i is estimated; The (7+) is 
an azimuthal average of the tangential component of the 
gravitational shear, which we take (7+(^)) ~ (.9+(^)) in 
the weak lensing limit. Then the projected mass of the 
cluster inside the projected radius 6 is given as 

Mc(<0) = ^02ScrCc(e; ^inn 7 ^out 

). (16) 

Mq{< 6) is regarded as a lower bound to the true enclosed 
mass M(< 0), because of the subtraction of the mean k, 
within the outer region, 0inn < 9 < Oout- Errors on Cc are 
calculated by propagating the rms errors (t+ (0„) (equation 
[14]) for the tangential shear measurement. 

The lensing properties such as k and ja depend on the 
source redshift and the background cosmology through 
the critical lensing surface mass density, Scr, defined in 
equation (4). For all of the cluster targets, we assume the 
mean redshift of the red background galaxies is (zs) = 1 
(e.g., Broadhurst et al. 2005). Since our target clusters 
are located at low redshifts {zd ^ 0.2 except Zd — 0.279 
for A1758) Ecr depends weakly on the source redshift as 
shown in Figure 2, so that a precise knowledge of the 
redshift distribution of background galaxies is not crucial 
(see, e.g., Bartelmann & Schneider 2001). 

In order to quantify and constrain the cluster mass 
properties we fit the following two different mass models 
to (^c-statistic measurements: singular isothermal sphere 
(SIS) and universal density profile proposed by Navarro, 
Frenk, & White (1996, hereafter NFW). The details of the 
halo models are given in Appendix. We parametrize our 
halo models using the cluster virial properties, such as the 
cluster virial mass, Mvir, the cluster virial epoch, Zvir, and 
the cluster virial radius, r^iT'. 

An 

Mvir = Yp(zvh.)Avh.r:i,, (17) 

where Avir is the overdensity with respect to the mean 
cosmic density /o(zvir) at the cluster virial redshift, pre- 
dicted by the dissipationless spherical tophat collapse 
model (Peeebles 1980; Eke, Cole, & Frenk 1996; Bullock 
et al. 2001). We assume the cluster redshift Zd is equal to 
the cluster virial redshift Zvh-- 

The SIS model has a one-parameter functional form de- 
scribed by the one-dimensional velocity dispersion of 
the cluster, which is related with the cluster virial redshift 
Zvir and the cluster virial mass Afvir as 

(Ti, (Afvir, ^^vir) = ^^-vir-ffo V^mO A^ir ( 1 + 2:vir)^ (18) 
1/3 

Thus. (Tt, cx Afy/. . Alternatively we can introduce the 
virial temperature of an SIS halo: 
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firupal = ksTsis (19) 

where /i = 0.62 is the mean molecular weight and nip is 
the proton mass. 

On the other hand, the NFW model has a two- 
parameter functional form, and we take the virial mass, 
-Mvir, and the concentration parameter, Cvir ^'"vir/rs, with 
Ts being the inner characteristic radius of the NFW pro- 
file. 

3. DISTRIBUTIONS OF MASS AND 
BARYONS IN MERGING CLUSTERS 

Here we present and compare in Figures 5-11 the re- 
sulting two dimensional maps of the convergence k, opti- 
cal luminosity density, and X-ray emission for our sample 
of merging clusters. Table 5 lists the properties of our 
background galaxy samples and the parameters relevant 
for weak lensing mass reconstructions. The cluster targets 
are selected to be on the various merging stages based on 
previous detailed X-ray studies, for which the physical re- 
lationship between mass and baryons though the merging 
processes is yet unknown. X-ray properties of the target 
clusters are summarized in Table 2. A1750 and A1758 
are binary clusters which are presumably in a pre-merger 
phase. A2034 and A2142 are cold front clusters. A754, 
A1914 and A520 are on-going mergers, which have irreg- 
ular temperature distributions and radio halos (Govoni ct 
al. 2004). 

3.1. A754 

A754 is an on-going merger at a redshift of 2 = 0.054. 
For such a low redshift cluster (2 < 0.1), the expected lens- 
ing signal is very low (Figure 2) and a weak lensing anal- 
ysis is challenging. As seen from the Chandra X-ray con- 
tours in the top-right panel of Figure 5, there are two ma- 
jor gas components in the east and west sides of the data 
field. The west gas clump (XC in Figure 5) corresponds 
to the main cluster having a BCG. From this configura- 
tion, the east X-ray gas clump (XE in Figure 5) seems 
to have been running through the cluster center from the 
west. The XMM-Newton image with a wider ficld-of-view 
has shown a moving feature of the east clump consistent 
with the Chandra X-ray observation (Henry, Finoguenov, 
& Briel. 2004). The X-ray contours of the east gas clump 
are compressed and elongated towards the north-east di- 
rection, indicating that the east gas clump is currently 
moving towards the northeast. Its averaged temperature 
within 9 arcmin radius is 10.0 ± 0.3 keV at the 90% confi- 
dence level (Markevitch et al. 2003; see also Table 2). The 
X-ray temperature maps of A754 were derived by Henry & 
Briel (1995) with ROSAT data, Henriksen & Markevitch 
(1996) with ASCA data, Markevitch et al. (2003) with 
Chandra data, and Henry, Finoguenov & Briel (2004) 
with XMM-Newton data. Henry & Markevitch (1996) 
obtained an X-ray temperature of 8.5-9.0 keV from the 
ASCA data. Henry et al. (2004) derived an X-ray tem- 
perature within 12' to be 8.6 ±0.1 ±0.6 where the first 
and the second errors are statistical and systematic un- 
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certainties at Icr confidence, respectively. The ASCA and 
XMM-Newton temperatures of A754 are somewhat lower 
than the Chandra temperature. However, the ASCA and 
Chandra measurements covered larger volumes of this on- 
going merger. The Chandra temperature map shows a 
strongly irregular feature, where the cool gas is offset from 
the X-ray brightness peak towards the northeast by ~ 200 
kpc. 

Since this nearby cluster extends outside the field- 
of-view of Subaru/Suprime-Cam. we used a finite-field 
method for reconstructing the mass distribution of A754 
as shown in Figure 5. We also compared the k map de- 
rived using the Kaiser & Squires method with the finite- 
field based k map. The k maps derived with the different 
methods are qualitatively similar, and the main features 
and quantitative properties of the mass peaks are consis- 
tent with each other within la reconstruction error, which 
means that the mass reconstruction is insensitive to the 
boundary conditions. The Gaussian FWHM used for the 
mass reconstruction is lf67. The reconstructed k map 
shows two mass clumps in the west and east sides of the 
data field. The west clump is located in the X-ray clus- 
ter central region (Henry, Finoguenov & Briel 2004), and 
hence we refer to this as the central mass clump (C in 
Figure 5). We found moderate signal-to- noise ratios of 
5.8 and 5.1 in the k map for the central and east mass 
(E in Figure 5) clumps, respectively, thanks to the superb 
image quality of Subaru/Suprime-Cam. 

The peak locations in the luminosity map of cluster 
sequence galaxies are in good agreement with these mass 
substructures (bottom- left panel) . The central mass peak 
coincides with the BCG (top-left panel). The east mass 
concentration contains a luminous, large elliptical galaxy. 

The east mass clump is offset from the cast X-ray clump 
towards the southwest, which is opposite to the moving 
direction of this gas clump (bottom-right panel). This 
kinematic feature is different from the case of the bullet 
cluster, 1E0657-56 (Clowe, Gonzale & Markevitch 2004). 
The cool gas region (see Markevitch et al. 2003) is located 
outside the region enclosed by the la mass contour around 
the east mass clump. 

A possible scenario explaining the observed merger ge- 
ometry of the east mass and X-ray clumps, which can- 
not be simply explained by the ram-pressure stripping, 
is the following: the east mass substructure just reaches 
its apocenter of the merger orbit and falls back towards 
the center for its second impact. In the rest frame of the 
mass clump, the X-ray core, which is initially bounded in 
the mass clump, feels the force in the opposite direction 
to the acceleration (and moving) direction of the mass 
clump. This explains the observed configuration of the 
mass and gas clumps associated with the east substruc- 
ture. As a result, the gas clump escapes away from the 
potential well of the mass substructure and would cool 
adiabatically as it expands. 

3.2. A 1750 

A1750 is a binary cluster at 2 = 0.086, well studied by 
Einstein, ROSAT, ASCA and XMM-Newton X-ray satel- 
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lites (Forman et al. 1981; Novicki, Jones, & Donnelly 
1998; Donnelly et al. 2001; Belsole et al. 2004). The 
XMM-Newton observation (Belsole et al. 2004) has shown 
that the southern cluster, A1750C, has a higher X-ray 
luminosity than the northern cluster, A1750N. The pro- 
jected distance between the two X-ray peaks is about 900 
kpc. The X-ray temperatures of A1750C and A1750N are 
3.87 ±0.10 keV and 2.84 ±0.12 keV at 90% confidence, 
respectively, (Belsole et al. 2004). The projected sep- 
aration is therefore shorter than the sum the two virial 
radii predicted by the M — Tx relation. Hence, it is likely 
that the two clusters just started to interact with each 
other. Indeed, the temperature in the middle region of 
the two clusters is higher than the cluster temperatures: 
Tx = 5.121°;^^ keV. Belsole et al. (2004) estimated a 
Mach number of 1.64 for the merger shocks by applying 
the Rankine-Hugoniot jump conditions under the assump- 
tion that the pre-shock temperature is approximated by 
the averaged temperature in the central regions of A1750C 
and A1750N and that the post-shock temperature is given 
by the temperature in the middle region of the two clus- 
ters. Belsole et al. (2004) also concluded that A1750C 
is an unrelaxed cluster for the following four reasons: (1) 
the discontinuity of the gas density profile in the south- 
east region (region 3 in their notation), (2) the shift of 
X-ray centroid from the position of the BCG, (3) the lack 
of evidence for a cooling flow, and (4) excess entropy in 
the central region, compared with other relaxed clusters. 

No significant offset among distributions of the galaxy 
luminosity, mass and ICM is found in the binary cluster 
A1750, as shown in Figure 6. 

The X-ray surface brightness peak in A1750C coincide 
with the position of the BCG (top-right panel) but is offset 
from the centroid of the X-ray emission, which is consis- 
tent with the X-ray analysis by Belsole et al. (2004). The 
position of the main mass peak {7.7a) in A1750C coincides 
with the peak positions of the optical luminosity (bottom- 
left panel) and the X-ray surface brightness (bottom-right 
panel), within a smoothing scale {1'.25 FWHM). No signif- 
icant mass substructure is found in A1750C. The observed 
geometry of mass, galaxies, and X-ray emission in A1750C 
does not show a signature of recent strong mergers as seen 
in A754. However, we cannot entirely rule out the possi- 
bility of a minor merger of small mass clumps below the 
weak lensing sensitivity. The averaged X-ray temperature 
of A1750C is higher than the temperature predicted by 
the best-fitting SIS model to the tangential shear profile, 
as shown in Figure 13 (see also §4.3). This is consistent 
with the lack of cool gas and the excess entropy in the 
central region (Belsole et al. 2004). 

A1750N has two BCGs whose positions coincide with 
that of the X-ray peak of A1750N (top-right panel). The 
mass peak (5.9ct) of A1750N is slightly, but not sig- 
nificantly, offset from the X-ray and galaxy-luminosity 
peaks, and the offset is within the smoothing scale of 
FWHM = 1.'25. The temperature Tsis estimated from 
weak lensing coincides with the X-ray temperature Tx 
for A1750N (see §4.3). 

As seen from the bottom-right panel of Figure 6, the 



reconstructed mass distribution traces an extended X- 
ray substructure (XM in Figure 6) between A1750N and 
A1750C, which is around (RA,DEC)= (202.77°, -1.765°). 
Overall, the ICM and galaxy distributions in the binary 
cluster A1750 and its components, A1750C and A1750N, 
are in good agreement with the mass distribution. 

3.3. A 1758 

A1758 is a binary cluster at z = 0.2790. The ROSAT 
observation by Rizza et al. (1998) showed that this cluster 
consists of two clusters (A1758N and A1758S) separated 
by ~ 8' on the sky, corresponding to the projected physi- 
cal separation of ~ 2Mpc. David & Kempner (2004) found 
from the XMM-Newton observation that the line-of-sight 
velocity difference between A1758N and A1758S is less 
than AV = 2100 kms"\ David & Kempner (2004) argued 
that their physical proximity to each other and small ve- 
locity difference are consistent with those of a gravitation- 
ally interacting system. The Chandra and XMM-Newton 
observations showed that there is no X-ray signature of 
the interaction between A1758N and A1758S (David & 
Kempner 2004) , contrary to the case of the binary cluster 
A1750 which has a high temperature region between the 
two components (see §3.2). This system is most likely in 
an earlier merger stage than A1750, and therefore serves 
as an ideal target for studying the cluster mass distri- 
bution in the initial stage of the merging process. As 
shown in Figure 7, the XMM-Newton X-ray image shows 
a bridge connecting between the two cluster components. 
On the other hand, the reconstructed mass map shows no 
significant mass structure corresponding to the bridge in 
the X-ray image. Note that the mass structure located 
on the X-ray bridge is a local minimum (M in Figure 7). 
It is also shown in Figure 7 that there is no significant 
offset between the mass structures and the X-ray/optical 
structures of A1758N and A1758S along the north-south 
direction connecting between A1758N and A1758S, which 
is consistent with results for A1750 (§3.2). 

The system components, A1758N and A1758S, are both 
undergoing mergers (David & Kempner 2004). David & 
Kempner (2004) found from XMM-Newton data average 
gas temperatures of A1758N and A1758S to be 8.2 ±0.4 
keV and 6.4+p 4 keV, respectively, at the 90% confidence 
level. Using Chandra data, David & Kempner (2004) 
obtained an average gas temperature of A1758N to be 
9.0~'^ j, keV at the 90% confidence level. For A1758N, the 
XMM-Newton and Chandra measurements of the average 
gas temperature are in good agreement with each other. 

A1758N has a complex X-ray morphology as seen in 
the top right panel of Figure 7: there are northwest 
and southeast X-ray subclumps in the cluster central re- 
gion. The Chandra X-ray image suggests that the north- 
west subclump (XC in Figure 7) is currently moving to- 
wards the north and that the southeast subclump (XSE 
in Figure 7) is currently moving towards the southeast 
(David & Kempner 2004). The optical luminosity dis- 
tribution (bottom-left panel) reveals two luminous sub- 
clumps of cluster-sequence galaxies in A1758N. The X- 
ray peak position of the northwest subclump coincides 



No. ] 

with that of the himinous BCG (top-right panel). The 
southeast gas subclump is, on the other hand, offset by 
w 290 kpc to the northwest of the other galaxy subclump 
(top-right panel). 

Figure 7 shows that mass and light are similarly dis- 
tributed in A1758N. The mass map shows double peaks: 
the first peak with a significance of 13.5ct (SE in Figure 7) 
corresponding to the southeast luminous galaxy clump, 
and the second peak with a significance of 11. Oct (C in 
Figure 7) corresponding to the northwest luminous galaxy 
clump (bottom- left panel). A visual inspection also re- 
veals several gravitational arc candidates around the two 
mass peaks in A1758N, supporting the bimodal mass dis- 
tribution in A1758N. Figure 15 displays zoom in views 
of tangential arc candidates in the i?c-band image. The 
tangential arc candidates A, B, and C are bluer than the 
cluster red sequence in g' — color. The tangential arc 
candidates D and E are associated with cluster galaxy con- 
centrations. This bimodal feature of A1758N in the weak 
lensing mass map was previously reported by Dahle et al. 
(2002) based on weak lensing data taken with ALFOSC 
on Nordic Optical Telescope. We note that Dahle et al. 
(2002) also found an arc-like image in A1758N. 

The northwest galaxy clump has a higher luminosity 
than the southeast one. The angular extent of the north- 
west galaxy clump is larger than that of the southeast one. 
These optical features would suggest that the northwest 
mass structure (C) is the primary component of A1758N 
and the southwest clump (SE) is the merging substruc- 
ture. TV-body simulations of ACDM models show that 
weak lensing measurements of cluster mass peaks have 
a large scatter (^ 100% for less massive systems with 
T\^200 ~ IO^^Mq/i"^) due to line-of-sight projection effects 
of intervening mass structures as well as due to intrinsic 
ellipticities of background galaxies (White, van Waerbeke 
& Mackey 2002; Wu et al. 2006). Therefore the rela- 
tive peak heights of the mass peaks SE and C (13.5cr and 
11.0(7 significance, respectively) could have been substan- 
tially affected by such effects. The southeast structure 
in the mass/galaxy map is located in front of the X-ray 
southeast subclump moving towards the southeast. This 
geometry regarding the southeast substructure is consis- 
tent with the results from the joint optical photometric, 
X-ray, and weak lensing analysis of the merging cluster 
1E0657-56 (Clowe, Gonzale, & Markevitch 2004). On the 
other hand, the northwest structure, which is likely to be 
the main component of A1758N, shows no significant off- 
set among the galaxy, ICM, and mass distributions, which 
is different from the results of 1E0657-56. 

A1758S has an elongated X-ray emitting core located 
close to the first peak in the galaxy luminosity distribu- 
tion (right panels). David & Kempner (2004) pointed 
out that the X-ray core has two bow-like shaped edges 
which are curved up toward north and south, respectively. 
They suggested that these features arc a signature of the 
merger along a line between the northwest and south- 
east sides with a small impact parameter. No optical 
counterpart has been detected in the Subaru image (left 
panels), whereas two luminous galaxy concentrations are 
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found along the northeast-southwest direction, which is 
perpendicular to the hypothetical merging direction. The 
second optical peak associated with A1758S is offset from 
the X-ray core to the southwest by ~ 480 kpc. The high 
density region with k > 9(t significance is elongated along 
the northeast-southwest direction, aligned with the two 
luminous galaxy concentrations. The mass peak asso- 
ciated with A1758S is located in the middle of the two 
optical clumps, and has a peak height of 9.9(7. The an- 
gular separation between the two luminous subclumps is 
comparable to the smoothing scale FWHM = 1.'25 for the 
weak lensing reconstruction. At this smoothing scale, we 
cannot therefore resolve possible mass substructures origi- 
nally associated with the optical subclumps. Nonetheless, 
the elongated mass and X-ray structures and the exis- 
tence of two luminous galaxy concentrations indicate that 
A1758S is not yet in dynamical equilibrium and is an on- 
going merger. The merger geometry of A1758S is puz- 
zling, because the curved-up directions of the X-ray core 
are almost perpendicular to both directions of the elonga- 
tions of the mass and light distributions. If the observed 
direction of the elongation in mass and light is regarded 
as the original collision axis, then our results would indi- 
cate that not all arc-like cores in the X-ray emission can 
be used to trace the path of the moving core. 

34. AI914 

The cluster A1914 is an on-going merger at z = 0.1712 
with averaged temperature of 10.9 ±0.7 keV at 90% con- 
fidence (Govoni et al. 2004). A1914 has two X-ray cores 
embedded in a round X-ray halo: As shown in the top- 
right panel of Figure 8 (see also Govoni et al. 2004) , the 
larger X-ray core (XE in Figure 8) is elongated in the east- 
west direction, and the smaller core (XW in Figure 8) is 
located to the west of the larger core. Based on Chandra 
X-ray data, Govoni et al. (2004) found a filamentary hot 
region along the northeast-southwest direction connecting 
between the two small cores, and argued that the filamen- 
tary hot region is a signature of the shock-heated gas in- 
duced by the elongated small X-ray core moving towards 
the west. 

We have obtained for A1914 a high resolution mass map 
with a Gaussian FWHM of 0'.75 as shown in Figure 8 (left 
panels). The cluster mass distribution is highly irregular, 
whereas the Chandra X-ray map shows a nearly circu- 
lar morphology except the central core region (top-right 
panel). We found seven mass concentrations above a de- 
tection threshold of Aa significance in the n map covering 
a 16' X 15'5 field around A1914. There are two promi- 
nent mass peaks with k 0.3: The first mass peak (CI in 
Figure 8) has a peak height of 9.1a, and its position is in 
good agreement with the second brightest cluster galaxy 
(BCG2 in Figure 8) which is close to the elongated X- 
ray core. The position of the second mass peak with a 
significance of 7.8a (C2 in Figure 8), on the other hand, 
coincides well with that of the brightest cluster galaxy 
(BCGl in Figure 8) sitting in the optical center. The 
observed geometry of the galaxy and mass distributions 
would suggest that the second mass peak (C2) associated 
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with the optical cluster center (BCGl) is the primary clus- 
ter center, and the first mass peak (CI) is the merging 
substructure. By visual inspection we found an arc-like 
image tangentially-orientcd with respect to BCG2/C1; as 
pointed out by Dahle et al. (2002) and Sand et al. (2005). 
The curvature of this tangential arc candidate is consistent 
with the existence of the high density peak (CI) around 
BCG2 (see Figure 16). 

However, it remains puzzling that the projected mass 
distribution in A1914 is so irregular contrary to the quasi 
circular X-ray distribution in the outer region. A detailed 
analysis will be presented in a forthcoming paper (Umetsu 
& Okabe in preparation). 

3.5. A 2034 

A2034 is a cold front cluster at z = 0.1130 well studied 
with ROSAT, ASCA and Chandra X-ray sateUites (David, 
Forman, & Jones 1999; White 2000; Kempner & Sarazin 
2003). Based on the Chandra observation, Kempner & 
Sarazin (2003) found a sharp discontinuity in the X-ray 
surface brightness ~ 3' to the north of the cluster center, 
similar to the surface brightness jump associated with cold 
fronts. The X-ray morphology suggests that this dense 
core is moving toward the north direction. Kempner & 
Sarazin (2003) showed that the temperature is fairly con- 
stant to a radius of at least about 5' from the center of the 
main cluster in A2034, and hence no significant tempera- 
ture jump, associated with the dense core, is found. They 
also argued that the excess emission from the southern 
part of the dense core is likely to be a background struc- 
ture. Kempner & Sarazin (2001) found evidence for the 
existence of a radio relic near the cold front. 

The resulting mass distribution in A2034 is highly irreg- 
ular and quite different from the X-ray surface brightness 
distribution of the XMM-Newton data as shown in Figure 
9. In the central region of this cluster (~ 18. '5 x 18. '5), 
there are six significant mass concentrations with peak 
heights greater than 4(t significance (left panels). The 
main mass peak with a significance of 6.4cr (C in Figure 
9) coincides with the first peak in the X-ray surface bright- 
ness distribution (see the bottom-right panel of Figure 9). 
There are two BCGs close to the first X-ray and mass 
peaks, corresponding to the main cluster of A2034. 

The second mass peak with a significance of 6.3(7 (S 
in Figure 9) is located to the south of the main mass 
peak, and coincides with the south excess X-ray emis- 
sion discovered by Kempner & Sarazin (2003). No spec- 
troscopic data are available in this region. The recon- 
structed mass map reveals another three mass structures 
(Wl, W2 and W3 in Figure 9) associated without X- 
ray counterparts in the west of the second mass peak. 
We confirmed that they are background structures at 
z ~ 0.116— 0.118 from SDSS spectroscopic data. The ve- 
locity difference between the main cluster and the west 
structures is cSz ~ 900— 1500 kms~^, and the comoving 
radial separation is C'(lOMpc) by neglecting their proper 
motions. Their apparent proximity suggests that the clus- 
ter A2034 and the west mass structures could be part of 
a filamentary structure. 



The primary result of our joint analysis on A2034 is 
the detection of the northern mass subclump with a sig- 
nificance of 5.4(7 (N in Figure 9) located ahead of the 
north cold front shown in the top-right panel. The sec- 
ond brightest member galaxy of A2034 is located close 
to the center of the north mass subclump (left panels). 
This merger geometry is consistent with the results of 
other cold front clusters, A2142 (§3.6), A1758N (§3.3), 
and 1E0657-56 (Clowe, Gonzale, & Markevitch 2004). 

3.6. A2142 

A2142 is a cold front cluster at z = 0.0909. Cold fronts 
were discovered for the first time in this cluster by the 
Chandra X-ray observation (Markevitch et al. 2000). 
There are two cold fronts in the central region of A2142, 
as indicated in the top-right panel of Figure 10: One is at 
the northwest edge (Cold Front NW) of a large core, and 
the other at the south edge (Cold Front S) of a small core 
in the X-ray surface brightness distribution. The small 
core is apparently contained within the large core in the 
sky plane. Those edges indicate that the large and small 
cores are moving toward the northwest and the south di- 
rections, respectively. The temperatures change abruptly 
across the cold fronts by a factor of two. There are two 
BCGs (BCGl and BCG2) in the cluster core region, as 
indicated in the top-right panel of Figure 10: BCGl is 
located inside the small core, and BCG2 is apparently 
located in the large core. The two BCGs have a large 
line-of-sight velocity difference of 1840 kms~^ (Oegerle, 
Hill, & Fitchett 1995). 

Okabe, Umetsu & Hattori (2008) have obtained a mass 
map with low angular resolution of FWHM = 3. '3 using 
a different set of Subaru weak lensing data taken under 
bad weather conditions. Okabe et al. (2008) found that 
the mass distribution is similar to the luminosity distri- 
bution of spectroscopically selected bright cluster mem- 
ber galaxies, but is quite different from the ICM dis- 
tribution. Okabe et al. (2008), based on the joint X- 
ray/optical/wcak-lensing analysis, found a marginal de- 
tection (2.1a) of a northwest mass substructure, which 
is located ahead of the northwest cold front. The pro- 
jected separation between the northwest cold front and 
mass substructure is ~ 540 kpc. 

In the present study we have derived a high resolution 
mass map with FWHM = IfOO, as shown in the left pan- 
els of Figure 10, using a new Subaru data set with long 
exposure. 

The number density of background galaxies usable for 
a weak lensing analysis is increased by a factor of ~ 3 as 
compared with the previous data set (Okabe et al. 2008). 

Using the new high-resolution mass map, this north- 
west mass structure (NW in the top-left panel of Figure 
10) has been detected at a significance of k = 3.0(7, which 
is located in the same region as the previous detection. 
This improved detection using an independent data set 
confirms our previous tentative detection, and supports 
the existence of the northwest mass substructure located 
ahead of the northwest cold front. In the bottom- left panel 
of Figure 10, we see a slight excess luminosity of clus- 
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ter sequence galaxies associated with the northwest mass 
substructure, as found earher by Okabe et al. (2008) us- 
ing the spectroscopicahy selected sample of bright cluster 
galaxies. The Chandra X-ray image shows no X-ray coun- 
terpart to the northwest mass substructure (right panels 
of Figure 10). XMM-Newton observations also show no 
substructure in the X-ray emission associated with the 
northwest mass substructure (Okabe et al. 2008). The 
merger geometry of A2142, that is the existence of mass 
substructures in front of cold fronts, is similar to other 
cold front clusters, such as A2034 (§3.5) and 1E0657-56 
(Clowe, Gonzale, & Markevitch 2004). 

We have also discovered a mass concentration associ- 
ated with a luminous galaxy concentration located in front 
of the south cold front. This south mass substructure (S 
in Figure 10) has a peak height of 4.3(T in the k map, and 
is located around (RA,DEC) = (239.54°, 27.18°). The line 
connecting between BCGl and this south substructure is 
roughly parallel to the moving direction of the south cold 
front. 

The new high resolution mass map has revealed more 
complex and detailed structures in A2142 than the low- 
resolution mass map by Okabe et al. (2008). The signif- 
icance of the first mass peak is 10.9cr in the k map, and 
the peak position is close to BCGl. We note the presence 
of two tangential arc candidates around BCGl as shown 
in panels A and B of Figure 17. 

The hne-of-sight velocity (LOSV) of BCGl is close to 
the central value (mode) of the LOSV distribution of clus- 
ter galaxies (Oegerle, Hill, & Fitchett 1995). The X-ray, 
optical, and mass peak positions of the main cluster are 
in good agreement with each other within FWHM = IfOO 
of the weak lensing reconstruction. On the other hand, no 
significant clumpy mass structure is seen around BCG2. 
Further we did not find any tangential arc candidate 
around BCG2. We show in panel C of Figure 17 a zoom 
in view of the 0f8 x 0.'8 region around BCG2. The pro- 
jected separation between BCGl and BCG2 is larger than 
the angular resolution of the mass map, FWHM = l.'OO. 
BCG2 has a large peculiar velocity, and is in the high- 
velocity tail of the LOSV distribution in the cluster rest- 
frame (Okabe et al. 2008). 

There are few member galaxies associated with BCG2 
in the LOSV distribution. The mass distribution around 
BCG2 is peculiar and different from all other BCGs in our 
target clusters, which are sitting at high density regions 
with large K-values. If BCG2 was originally associated 
with the northwest mass substructure, then the results 
may suggest the possibility that BCG and the dark-matter 
halo underwent different dynamical processes, which then 
resulted in a large offset between the BCG and the dark- 
matter halo after the merger. 

3.7. A520 

Markevitch et al. (2005) discovered with Chandra X- 
ray observations a bow-shaped shock close to a dense gas 
clump in the cluster A520 at z — 0.199. They derived a 
Mach number of M = 2.1^03 with the Rankinc-Hugoniot 
condition. The post and pre-shock temperatures are found 
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to be 11.5^11 keV and 4.8tj^J keV, respectively. The X- 
ray feature is very similar to the bullet cluster 1E0657-56 
(Markevitch et al. 2005), in that a shock detached from a 
dense core is prominently observed. Therefore this cluster 
is a good target for a weak lensing study of 'bullet '-type 
merging clusters. 

Two compact gas clumps (XSWl and XSW2 in Figure 
11) are located in the southwest of the cluster center. 
The gas clump XSW2 associated with the bow shock is 
close to but slightly offset from a luminous BCG. The 
extended tail gas from XSW2 is also seen in Figure 11. 
Markevitch et al. (2005) argued that XSW2 is moving 
along the northeast-southwest direction. Govoni et al. 
(2004) found that the temperature is higher along this di- 
rection based on Chandra data. On the other hand, the 
gas clump XSWl is elongated in the east-west direction, 
and its X-ray contours are compressed towards the north- 
west direction (see the right panels of Figure 11), suggest- 
ing that XSWl is currently moving towards the northwest 
direction. The presence of two gas clumps might suggest 
that A520 is a three-component merger. 

Our weak lensing analysis of A520 is based on archival 
i'-band data retrieved from SMOKA. We have three i' 
imaging data sets used for our weak lensing analysis as 
described in §2.1. Here we shall focus on the analysis 
results using the deepest data set (see Table 1), which is 
based on the same data used by Mahdavi et al. (2007). 
A comparison of mass maps with the three different data 
sets is devoted in Appendix 2. 

We show in the left panels of Figure 11 contour plots 
of the reconstructed mass distribution in A520. The 
smoothing scale for the weak lensing mass reconstruction 
is FWHM = 1.'25. The first mass peak (CI in Figure 11) 
is located in the southeast of the optical cluster center, 
and the detection signal-to-noise ratio is 6.1a. The mass 
peak CI is located ~ 0'.5 east of the central X-ray surface 
brightness peak (XC in Figure 11). Recently Mahdavi et 
al. (2007) reported the presence of a massive dark core 
(the mass peak 3 of Mahdavi et al. 2007) located ~ 0.'7 
south of the central X-ray peak, XC; the offset between 
the dark core and XC is smaller than the angular resolu- 
tion of Mahdavi et al.'s mass map (FWHM = 1'). The first 
peak CI in our mass map and the dark core of Mahdavi 
et al. (2007) coincide with XC within Gaussian FWHMs 
of our and their mass maps, and the projected distance 
between CI and the dark core is about I'.O. In the cen- 
tral region of the n map, there are two more significant 
local maxima (C2 and C3 in Figure 11): C2 with a signif- 
icance of 5.8(7 is located in the north of CI, and C3 with 
a significance of 5.4(t in the east of CI. The locations of 
C2 and C3 agree well with those of the main and a lo- 
cal maxima of the optical luminosity density distribution 
(see the bottom-left panel of Figure 11), respectively. In 
the central high-density region of A520, the shape of mass 
contours with k ^ 4cr significance is similar to that of the 
cluster galaxy luminosity distribution. 

We note that although Mahdavi et al. (2007) used the 
same Subaru/Suprime-Cam i' data as in the present pa- 
per, their analysis method is different from ours in the fol- 
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lowing way: Our weak lensing analysis is based on three 
sets of stacked Subaru/Suprime i' data (see Appendix 2), 
while Mahdavi et al. (2007) performed shape measure- 
ments separately for each of seven Subaru exposures and 
combined individual weak lensing catalogs to improve the 
accuracy for the shear estimates. We have confirmed the 
existence of the dark core (CI) reported by Mahdavi et al. 
(2007) in two of the imaging data sets; however, the peak 
CI is not particularly pronounced in the mass map based 
solely on i' images taken without AG. We have detected 
C3 at high significance levels ( ^ 4.5(t) from all of the three 
imaging data sets (see Appendix 2), while Mahdavi et al. 
(2007) found no significant mass structure around C3, 

The second mass peak (SW2 in Figure 11) with a signif- 
icance of 6. Oct is associated with the southwest substruc- 
ture. In the southwest region, mass and light are simi- 
larly distributed, and the peak positions in the mass and 
luminosity maps coincide well with each other (see the 
bottom- left panel of Figure 11). On the other hand, the 
X-ray surface brightness peak associated with the south- 
west X-ray core (XSW2 in Figure 11) is slightly offset from 
the corresponding galaxy-luminosity and mass peaks; the 
separation between the X-ray and mass peaks is smaller 
than the Gaussian smoothing scale of FWHM = 1.'25 used 
for the weak lensing mass reconstruction. The southwest 
compact gas, XSW2, seems to be within the potential 
well of the southwest mass substructure, SW2. We also 
find a local mass peak with a significance of 3.5ct (SWl 
in Figure 11) located ahead of the northwest of the gas 
clump XSWl. The mass peak SWl is associated with a 
concentration of cluster member galaxies (see the bottom- 
left panel of Figure 11). The projected separation between 
SWl and the northwest edge of XSWl is about 1.'6. 

We note that, in addition to these mass peaks, the re- 
constructed K map shows three more local maxima above 
a significance of 4ct, which are indicated as N, NEl, and 
NE2 in Figure 11: the mass clumps NE2 and N are associ- 
ated with slight concentrations of cluster member galaxies. 
On the other hand, there is no apparent concentration of 
galaxies having similar {V — i') colors around the location 
of NEl, although this mass structure is seen all of the 
three data sets (see Appendix 2). 

Finally, the results from our joint optical, X-ray, and 
weak-lensing analysis show that although the X-ray fea- 
tures of A520 are very similar to those of the bullet cluster 
1E0657-56, the relative positions of the galaxy. X-ray, and 
mass substructures are different, in that there is no appar- 
ent offset between the mass (SW2) and the gas (XSW2) 
clumps associated with the bow shock. 

4. GLOBAL CLUSTER PROPERTIES 

We present and compare global cluster properties, such 
as the total mass, optical luminosity, and X-ray tempera- 
ture, for our sample of merging clusters based on the joint 
optical-photometric/X-ray/weak-lensing analysis. 



4-.1. Cluster Mass 

We show in Figure 12 radial profiles of the projected 
mass Af^ (9) for our sample of seven merging clusters to- 
gether with the best-fitting NFW and SIS profiles. Note 
that the error bars are correlated; in the model fitting, 
we take into account the full error covariance matrix for 
the Mq measurements. The projected mass M/^iO) is 
measured using aperture densitometry, Cc, by M^(< 6) = 
TT{Dd9)^i:,,Cc{0;9inn,0out) (scc equation [15] and §2.3.4). 
Table 7 summarizes the best-fitting SIS and NFW model 
parameters derived for our sample of merging clusters. 
For each target and its cluster component (if exists) we 
chose the position of the BCG as the cluster center for 
measuring global cluster properties. The 0inn and 6'out 
define the inner and outer boundaries of the background 
annulus for the Cc-statistic measurement. The values of 
(^innj^out) ciTe choscu SO as to avoid cluster binary com- 
ponents and/or significant background substructures, and 
they are summarized in Table 8. The radial range used for 
the fit is different for each cluster: the upper limit, ^max, is 
set to the maximum radius which will avoid contamination 
by neighboring substructures and foreground/background 
mass structures; the lower limit, ^min, is set so as to avoid 
the strong-lensing regime (see Table 7). As we have seen 
in §3, the mass distribution in merging clusters is highly 
irregular and complex. Indeed, when fitting NFW/SIS 
models to the (^c-statistic measurements, we found for 
some clusters large x^-values relative to the degrees-of- 
freedom (dof). For all of the cluster targets, the virial 
mass Mvir of the best-fitting NFW model is consistent 
with that of the SIS model (sec §2.3.4) within Ict statisti- 
cal error. On the other hand, the NFW concentration pa- 
rameter Cvir was poorly constrained, except for a marginal 
detection of A2142 (w 2.5ct), because of contamination 
by neighboring substructures and clusters in projection 
space. 

As defined by equation (15), M^(< 9) yields a lower 
bound to the true enclosed mass M{< 9) due to the 
subtraction of the mean k over the background annulus, 
'*(^inn < 9 < 0out). In Table 8 we summarize the results 
of weak-lensing total mass estimates, A/^(< 9m), for our 
sample of merging clusters. The aperture radius 9m is set 
to the maximum radius available for each target cluster, 

i.e. 9m ~ ^'max- 

4-.2. Mass-to-light Ratio 

We measure for individual targets the cluster luminosity 
from color-selected member galaxies (see §2.2). The clus- 
ter luminosity is measured in i?c-band for all clusters, ex- 
cept in j'-band for A520. Here the aperture and the back- 
ground annulus for each target are taken to be the same as 
in aperture densitometry. Table 8 lists the mass-to-light 
(M/L) ratio derived for our sample of merging clusters, 
where the errors include a systematic uncertainty in the 
background luminosity density estimation as well as Ict 
random uncertainties in the cluster luminosity and mass 
measurements; the lower limit on the cluster M/L ra- 
tio is based on the cluster luminosity without background 
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subtraction and la random uncertainties for the M/L ra- 
tio. The clusters A754 and A1750 have considerably high 
M/L ratios of M/L > 600/i(M/L)o, which is simply due 
to the fact that cluster galaxies are still dominating in 
the background annular region and hence the background 
luminosity density is overestimated. 

Girardi et al. (2000) studied optical properties of a 
sample of 105 nearby Abell-ACO clusters (z < 0.15) for 
which virial mass estimates are available, and found typ- 
ical values for M/Lb^ ratios to be 230 — 250ft-(Af/iB^. )©. 
Sanderson & Ponman (2003) conducted a joint optical/X- 
ray study of the mass composite of a sample of 66 relaxed 
groups/clusters. Based on the X-ray mass measurements, 
they obtained a logarithmic mean value of M/Lb ~ 
350h{M/L)Q, which is scaled to M/Lr^ sa 190h{M/LRjQ 
assuming (L/Lq) b^ = 0. 547 {L/Lq)^^ from a typical color 
{Bj — R) = 1.702 of early- type galaxies (equation [4] of 
Sanderson & Ponman 2003). 

Clowe et al. (2006) performed a weak lensing analy- 
sis of 20 high-z clusters (z > 0.4) selected from the ESQ 
Distant Cluster Survey (EDisCS) using deep 3-color op- 
tical images taken with the VLT/F0RS2. They derived 
for an individual cluster M/L ratios, in / and B pass- 
bands, within 500 kpc from photo-z selected cluster galax- 
ies. They found for their optically-selected EDisCS clus- 
ters that the clusters tend to have a lower AI/L ratio at 
higher redshifts, but found no change in their M/L ra- 
tios with cluster mass. Clowe et al. (2006) also showed 
that the clusters having additional structures in projec- 
tion space have a higher measured M/L ratio, which is 
probably due to the projection effects. 

4-3. X-ray vs. SIS Temperatures 

Wc compare for individual cluster targets the X-ray 
temperature, Tx, and the best-fitting SIS temperature, 
Tsis, obtained from the aperture densitometry measure- 
ments. In Table 8 we summarize the best-fitting SIS ve- 
locity dispersion, cr^, from the weak lensing analysis. The 
X-ray temperatures of the clusters are taken from the lit- 
erature (see Table 2). We note there is a good agreement 
between the Chandra and XMM-Newton temperatures of 
A1758N (David & Kempner 2004), although Chandra and 
XMM-Newton have different instrument responses as well 
as sky and instrumental background levels. The X-ray 
temperature measurement is performed within apertures 
of 6* ~ 10', while the weak lensing measurement of aper- 
ture densitometry within apertures of ~ 5' — 15' (see 
Table 8), depending on the availability of the clean back- 
ground aperture without any significant mass structures. 
We emphasize that the SIS temperature Tsis based on 
gravitational lensing does not depend on the dynamical 
and physical state of the cluster system. 

Figure 13 compares the X-ray temperature and the SIS 
temperature for our sample of merging clusters. It is 
shown in Figure 13 that there is an overall trend that the 
ICM temperatures from the X-ray observations are higher 
than the SIS temperature from the weak lensing analysis. 
Figure 14 shows explicitly the X-ray to SIS temperature 
ratio, Tx/Tsis, as a function of Tsis. The sample mean 



for Tx/Tsis is estimated as (Tx/Tsis) =2.05 ±0.11. 

For A1750N, where mass, light, and X-ray emission 
are similarly distributed (see §3.2), the temperature ra- 
tio Tx/Tgis is consistent with unity within the la error 
bar. The ICM temperature of A1750C, which shows an 
excess entropy in the central region, is higher than Tsis 
(§3.2), but is consistent with Tsis within 2a: Tx/Tsis = 
1.72 ±0.53 (la). The averaged X-ray temperature of 
A520, Tx = 7.1 ±0.54 keV (la), is slightly higher than 
the SIS temperature, Tsis — 5.9 ± 1.1 keV (la), but is 
marginally consistent within la with Jsis^ Tx/Tsis = 
1.18 ±0.25 (la); However, we note that the post-shock 
temperature of A520, Tx = ll-5i|:i keV (90%CL), is sig- 
nificantly higher than Tsis? while its pre-shock tempera- 
ture, Tx=4:.8tli kcV (90%CL), is marginally consistent 
with Tsis. 

For the rest of the sample (A754, A1758N, A1758S, 
A1914, A2034, A2142) consisting of on-going mergers and 
cold-front clusters, the ratios Tx/Tsis are in the range of 
1.4-3.1, with the sample mean of (T^/Tsis) = 2.41±0.16. 
Hence, presumably the most promising scenario is that the 
ICM temperature is increased by heating processes trig- 
gered by cluster mergers. On the other hand, the temper- 
ature ratios for the cold front clusters, A2034 and A2142, 
are in a moderate range of 1.94 ±0.38 and 1.41 ±0.19, 
respectively. 

Although the averaged temperatures Tx of the cold- 
front clusters were measured including the cold, dense 
cores, the X-ray temperatures outside the dense cores 
agree with the global temperatures (including the cold 
cores) within errors (Markevitch et al. 2000; Kempner 
& Sarazin 2003). Therefore, the observed variations and 
trends in the Tx /Tsis ratios would be naturally explained 
by the different energy scales injected into the ICM or 
different merger states of the cluster mergers. Figure 14 
shows no clear correlation between the temperature ratio 
Tx/Tsis and the SIS temperature Tsis. This might sug- 
gest that the enhancement of the ICM temperature does 
not depend strongly on the mass of the main cluster but 
on merging stages and conditions. 

5. DISCUSSION AND CONCLUSIONS 

We presented and compared weak lensing mass, optical 
light, and X-ray emission maps for a sample of seven merg- 
ing clusters of galaxies. We have selected for the present 
study seven nearby Abell clusters (0.0542 < z < 0.279) of 
different merging stages and properties, and conducted 
systematic, deep imaging observations of the seven tar- 
get clusters with Suprime-Cam on Subaru telescope. Our 
seven target clusters, representing various merging stages 
and conditions, allow us to investigate in details the phys- 
ical interplay between dark matter, ICM, and galaxies as- 
sociated with hierarchical structure formation. The clus- 
ters A1750 and A1758 are binary systems, each of which 
consists of two cluster-sized components and has not 
yet experienced the first impact. A520, A754, A1758N, 
A1758S, and A1914, on the other hand, arc classified as 
on-going merging clusters. A2034 and A2142 are "cold 
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front" clusters. 

In the binary cluster A1750, the projected mass, op- 
tical light, and X-ray distributions all revealed clear bi- 
nary structures associated with the two cluster compo- 
nents (A1750N and A1750C), and are overall similar and 
regular without significant substructures (§3.2). A1758 
also exhibits binary structures associated with the two sys- 
tem components (A1758N and A1758S) in the projected 
mass, optical light, and X-ray distributions (§3.3). Along 
the line between A1758S and A1758N, which is parallel to 
the dynamical motion of the binary components, no signif- 
icant offset is found between the mass, X-ray, and optical 
peaks in the binary components. On the other hand, the 
ICM distribution in individual components, each of which 
is an undergoing merger, appears to be different from the 
distributions of mass (dark matter) and light (galaxies). 
Our results on the binary clusters indicate that, in an early 
merging phase, the distributions of mass and baryons are 
overall similar and the main mass peak of the cluster coin- 
cides well with the corresponding X-ray peak and galaxy 
concentration. 

The mass distributions of on-going merging clusters, 
A520, A754, A1758N, A1758S, and A1914, are found to 
be highly irregular. Overall, the mass distribution ap- 
pears to be similar to the galaxy luminosity distribution, 
whereas their distributions are totally different from the 
ICM distribution. Discrepancies between mass and ICM 
distributions are different in a variety of ways. Our re- 
sults on the on-going mergers can be formally classified 
into three types depending on the moving direction and 
relative positions of mass and X-ray substructures. Here 
wc speculate on the moving direction of cool, dense cores 
from the direction of compression and elongation of X-ray 
cores and trails, and summarize the results as follows: 

(1) No significant offset between X-ray and mass peaks 

The mass peak position coincides with the peak position 
of X-ray surface brightness within Gaussian smoothing 
FWHM of the weak lensing mass reconstruction. This 
apparent feature would imply that these dense cores are 
trapped in the gravitational potential well of the mass 
concentration. This is found in the southwest clump (SW2 
and XSW2) of A520, the northwest clump (C and XC) of 
A1758N, and the first mass peak and the elongated gas 
core (CI and XE) of A1914. The position of the BCG in 
A1758N coincides with the peak position of the dense core 
in the ICM (XC). For the case of A520 and A1914, the 
associated luminous galaxy is slightly offset perpendicular 
to the moving direction of the gas core. 

(2) Mass clump in front of the dense core 

The second type is the case where the mass clump is 
located ahead of the moving gas core, as seen in the south- 
west structure (SE and XSE) in A1758N. The bullet clus- 
ter 1E0657-56 is categorized into this type. We note that 
this configuration is also found in our cold front clusters, 
A2034 (N and the cold front) and A2142 (NW and S), as 
we will discuss below. The following scenario could ex- 
plain this type of merger configurations: a dense gas core, 



which is originally bound in a merging substructure, is 
stripped away from the mass and galaxy components of 
the substructure by the ram pressure of the ICM. 

(3) Mass clump behind the dense core 

The third type is where the mass clump is located be- 
hind the dense core, which is opposite to the second case. 
This is found in the east substructure (E and XE) of A754. 
A possible explanation is that this cluster is in a merging 
phase right after reaching its apocenter of the merger or- 
bit (see §3.1), although no direct evidence from currently 
available data. 

Cluster mergers depend on a certain set of merging pa- 
rameters, such as the mass ratio of the main to sub clus- 
ter, initial velocity and angular momentum of the merg- 
ing substructure, and the shape of the gravitational po- 
tential of the main cluster. Furthermore, clusters are in 
various phases of the merging process. Understanding 
and constraining the dynamical process that causes offsets 
between the mass/ICM/galaxy components therefore re- 
quire a more detailed, quantitative study of these merged 
substructures, which will be presented elsewhere. For 
example, line-of-sight kinematic information of member 
galaxies from spectroscopic observations, measurements 
of the density, temperature, and entropy of the pre/post- 
shock regions will be valuable. 

The mass maps of the cold front clusters A2034 and 
A2142 reveal irregular mass structures, which are quite 
different from the ICM distribution. Mass structures with 
low gas mass fractions have been detected in front of three 
cold fronts, which is found earlier in the bullet cluster 
1E0657-56 (Clowe et al. 2004). Such merger configura- 
tions are classified as the second type of on-going mergers 
as discussed above. To date four cold fronts have been 
studied via joint X-ray/weak- lensing analysis including 
1E0657-56, and all of the cold fronts show such a config- 
uration of mass/cold-front substructures. This might be 
thus a common feature of cold front clusters associated 
with the formation of cold fronts. Current results from 
the joint analysis of cold front clusters suggest a possible 
scenario for the formation of cold fronts that a substruc- 
ture falls into and passes through the primary cluster halo, 
while the hot gas originally bound in the substructure is 
stripped away by ram pressure stripping (see Clowe et al. 
2004). A more detailed, statistical study of the relation- 
ship between cold fronts and associated mass structures 
is required to draw more definitive conclusions on the for- 
mation of cold fronts. 

The observed positions of dense cores in A2034 and 
A2142 are apparently located in the primary cluster halo 
in projection space. This is different from the case of the 
bullet cluster 1E0657-56, which shows a significant offset 
between the mass and X-ray halos of the primary clus- 
ter as well as between those of the subcluster (Clowe et 
al. 2004). However, since no line-of-sight information of 
substructures is available from the current data sets, the 
merger geometry in three-dimensional space has not yet 
been fully constrained. Such line-of-sight information of 
the primary and sub clusters will be necessary to place 
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quantitative constraints on the evolution of the gas tem- 
perature and pressure of the ICM substructure against the 
gravitational potential of the primary cluster. 

Our quantitative comparison of projected distributions 
of mass, galaxies, and ICM in merging clusters shows 
that the dark matter and ICM components exhibit differ- 
ent behaviors during the cluster merger process. Similar 
features have been reported based on cosmological A^- 
body/hydrodynamical simulations. Tormen, Moscardini, 
& Yoshida (2004) have shown by their numerical simu- 
lations that dark matter halos of merged substructures 
freely move and oscillate in the primary cluster halo, while 
gas halos initially bound in the substructure are stripped 
by various physical processes, such as ram pressure strip- 
ping and dynamical friction. These behaviors depend on 
the mass ratio between the primary and sub cluster. It 
is easier for dark matter of small substructures to survive 
longer, while for the ICM of small substructures to be 
destroyed or digested by the primary cluster. They also 
showed that the separation between the dark matter and 
gas centers for a merged substructure, whose mass is more 
than 0.01 times of the primary cluster mass, is very small 
in the early phase of ^ 1 Gyr. These simulation results 
are consistent with our results of the binary cluster A1750 
(§3.2). Takizawa (2006) studied the X-ray and mass distri- 
butions of 1E0657-56 based on A^-body/hydrodynamical 
simulations, and reproduced a clear offset between X-ray 
and mass peak as found from the joint X-ray/ weak- lensing 
analysis by Clowe et al. (2004). The mass ratio of the sub 
to main halo is assumed to be 1/16, which is based on 
the weak lensing mass estimate by Clowe et al. (2004). 
Takizawa (2006) also derived an analytic expression for 
the ram pressure stripping conditions of the substruc- 
ture in mergers of two NFW dark matter halos, which 
is useful to understand the merging conditions of the first 
and second types defined above. Rickcr & Sarazin (2001) 
and Mathis, Lavaux, Diego, & Silk (2005) studied clus- 
ter mergers of nearly-equal-mass with numerical simula- 
tions and showed that the hot gas of a merged substruc- 
ture escapes from the local gravitational potential well 
when it reaches the apocenter in merger orbits, because 
the inward motion of the gas is delayed compared with 
the associated dark matter clump. Their simulation re- 
sults would be a clue to understand the third type of 
on-going mergers as found in A754. Nagai & Kravtsov 
(2003) searched for cold fronts and their counterparts in 
their high-resolution cluster simulations, and found cold 
fronts appear to be common in major and minor clus- 
ter mergers in hierarchical models. Their results indicate 
cold fronts are non-equilibrium transient phenomena dur- 
ing cluster mergers, although no significant offset between 
dark matter and gas distributions is shown in their figure. 
Ascasibar & Markevitch (2006) used high-resolution sim- 
ulations of cluster mergers and showed that cold fronts 
in relaxed clusters are due to sloshing of the cool gas in 
the gravitational potential of the primary cluster, and can 
be easily triggered by minor mergers and persist for gi- 
gayears. Based on three-dimensional hydrodynamic sim- 
ulations Takizawa (2005) investigated the dynamical evo- 
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lution of the ICM driven by a radially-moving substruc- 
ture and its observational implications, and found that 
the subcluster's cold gas is pushed out of its potential 
well around the turn around, resulting in that the cold 
gas clump appears to be in front of the main X-ray peak 
of the subcluster (see Figure 6 of Takizawa 2005). 

A comparison of the ICM and SIS temperatures of 
merging clusters from X-ray and weak-lensing analyses, 
respectively, provides us with interesting pieces of evolu- 
tionary information of the ICM during the cluster merger. 
Our results show that the ICM temperature in on-going 
and cold-front clusters is significantly higher than the SIS 
temperature that is a measure of the virial temperature of 
the cluster (Figure 13). This would suggest that the ICM 
temperature is increased by the cluster merger. Similar 
results were found from numerical simulations. Randall, 
Sarazin, & Kicker (2002) and Rowley, Thomas, & Kay 
(2004) studied the effects of merger boosts on the ICM 
properties and showed that cluster mergers can boost for 
a duration of the sound crossing time the X-ray luminos- 
ity and temperature of the merged cluster above the virial 
equilibrium values for the merged system. Our results 
show no significant correlation between the SIS tempera- 
ture, or the cluster mass, and the X-ray to SIS tempera- 
ture ratio, Ta'/^sis; this might indicate that the energy 
release from cluster mergers is determined not only by the 
cluster mass but also by other merging parameters. Rickcr 
& Sarazin (2001) studied offset mergers between clusters 
based on iV-body/hydrodynamic simulations as a func- 
tion of impact parameter and mass ratio of the colliding 
clusters, and show that the variations in global properties 
of clusters, namely total luminosity and average temper- 
ature, indeed depend on the impact parameter and the 
mass ratio, as well as the observed epoch of mergers. 

Based on visual inspection, we have identified gravi- 
tational arc candidates in three of our cluster targets: 
A1758, A1914, and A2142 (see Figures 15, 16, and 17). 
A further examination and confirmation of the lensing 
hypothesis of these gravitational arc candidates will re- 
quire follow up spectroscopy and higher-resolution imag- 
ing, which will provide strong lensing constraints needed 
for a detailed modeling of the observed arc-cluster sys- 
tems. 

In summary we compared projected distributions of 
mass and baryons in a sample of seven merging clus- 
ters based on the joint optical-photometric/X-ray/weak- 
Icnsing analysis. The global cluster parameters, such as 
the cluster mass, the cluster mass-to-light ratio, and the 
ICM temperature, are derived for individual targets of our 
cluster sample. 

A joint analysis of optical photometric. X-ray, and weak 
lensing data provides us with a comprehensive picture of 
cluster mergers, which will help to better understand the 
complex physical interplay between dark matter, ICM, 
and galaxies in the cluster formation process. For merg- 
ing clusters an X-ray analysis alone cannot constrain the 
cluster mass distribution since the ICM is not in hydro- 
static equilibrium. We have demonstrated that weak lens- 
ing is indeed a unique and powerful method to map the 
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mass distribution in merging clusters that are in the pro- 
cess of active evolution and formation, since it does not 
require any assumption of the physical/dynamical state of 
the system. Furthermore, thanks to the Subaru/Suprimc- 
Cam with excellent image quality and wide-field imaging 
capability, we were able to achieve high-resolution and 
wide- field imaging in weak lensing mass reconstructions. 
In particular, it is worth noting that the Subaru/Suprime- 
Cam makes it possible to reveal a mass distribution of a 
very nearby cluster at z ^ 0.05 (A754), where the lensing 
signal is very low due to its proximity to the observer. A 
detailed study of individual targets using the joint data 
sets is the next step for understanding the merger dy- 
namics of clusters and associated physical processes of 
the ICM. To date galaxy clusters have been studied ex- 
tensively and systematically in X-ray and optical photo- 
metric observations. On the other hand, the weak-lensing 
database for clusters is not sufficiently large as compared 
with X-ray and optical, although much progress has been 
made in the recent years (e.g., Dahlc et al. 2002; Cypriano 
et al. 2004; Bardeau et al. 2005, 2007; Clowe et al. 2006; 
Hoekstra 2007). 

Therefore a systematic weak-lensing study of a larger 
sample of clusters needs to be conducted for further un- 
derstanding of clusters and their formation/evolution pro- 
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Appendix 1. Lensing Mass Models 



and the virial mass Afvir of the SIS halo as 



The SIS density profile is given by 



(Al) 



where Ci, is the one-dimensional velocity dispersion of the 
SIS halo. The ay is related with the virialization epoch 



vir : ■^vir / 



' vir 



1 
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47r _ 

— p(Zvir)Av 



(A2) 
(A3) 



where rvir = 7'vir(A'/vir, Zvir) is the virial radius, and Avir 
is the mean overdensity with respect to the mean cos- 
mic density p(zvir) = ^miz^ii) Pcrizy-n-) at the virialization 
epoch, predicted by the dissipationless spherical tophat 
collapse model (Peeebles 1980; Eke, Cole, & Frenk 1996; 
Bullock et al. 2001). We assume the cluster redshift z is 
equal to the cluster virial redshift Zvir. We use the fol- 
lowing fitting formula in a flat 3-space with cosmological 
constant (see Oguri, Taruya, Suto 2001): 

(A4) 



Avir = 187r"(l + 0.40930;;;:^^"), 
where o^vir = l/ilm{zvii) — 1. 

1 /3 

We thus have cr„ oc M^,/. . In the case of an SIS lens, the 
lensing convergence is obtained as 

«sis(e) = ^, (A5) 

where Oe = 47r((T„/c)^(Dds/-Ds) is the Einstein radius. 
Then the averaged convergence within 9 is 



2Aesis(^). 



(A6) 



The NEW universal density profile has a two-parameter 
functional form as 

Ps 



PNFw(»') 



(A7) 



(r/r,)(l + r/r,)2 

where ps is a characteristic inner density, and rg is a char- 
acteristic inner radius. In stead of using Ts, we introduce 
the concentration parameter, Cvir = rvir/i's ■ The inner den- 
sity Ps can be expressed in terms of other virial properties 
of the NEW halo: 



- /O(Zvir) 



Av 



3 ln(l + Cvir)-Cvir/(l + 



(A8) 



Per is the critical density of the universe at the cluster 
redshift Zd, and rg is the NEW scaling radius. Hence, 
for a given cosmological model and a halo virial redshift 
(zvir), we can specify the NEW model with the halo virial 
mass Mvir and the halo concentration parameter Cvir. 

Eor an NEW profile, it is useful to decompose the con- 
vergence k{6) and the averaged convergence k{< 9) as 



%Fw(< x) ^ -g{x), 



(A9) 
(AlO) 



where h = ApsTg /'Sict is the dimensionless scaling conver- 
gence, x~9/ {rs/Dd) is the dimensionless angular radius, 
and f{x) and g{x) are dimensionless functions. We have 
analytic expressions for J{x) and g{x) as (Bartelmann 
1996): 



No. 



Subaru Weak Lensing Study of Seven Merging Clusters 



17 



1 



-1 



^^arctanh^i^ 
(x = l) 



(x<l) 

(All) 
(x>l). 

(x<l) 
(x = 1) (A12) 

Finally, we use the following identity to calculate the 
theoretical ^c-statistic for a given set of the aperture pa- 
rameter (6'in„,6lout): 




=arctan 



C,(0) = s(<0)-k( 



< 6) < 



t) 



k{< 9) - k{< Oinn) 
1 



(A13) 

(A14) 



'1 ^fl ^(S(<0in„)-K(<eout)).(A15) 

In practice, we use a discretized estimator for equation 
(15). We note that, unlike tangential shear measurements, 
measurement errors for Cc(^) arc correlated between dif- 
ferent annular bins. 

Appendix 2. Comparison of Weak Lensing Mass 
Reconstructions between Three DiflFerent Data 
Sets for A520 

We retrieved from the Subaru Archival Website, 
SMOKA, a total of seven I'-band images of A520 taken 
under good seeing conditions (~ 0."6). The archival i' data 
of A520 were taken with and without AG (acquisition and 
guide) probe on 17th November 2001 and 19th October 
2001, respectively. The use of AG probe is required for 
accurate telescope pointing, which is crucial for weak lens- 
ing shape measurements of faint background galaxies. We 
performed a weak lensing analysis separately on the fol- 
lowing three imaging data sets: (A) all of seven i' im- 
ages taken both with and without guide probe (7 x 240s 
exposure) used in the weak lensing analysis of Mahdavi 
et al. (2007), (B) four i' images taken without guide 
probe (4 x 240s exposures), and (C) three i' images taken 
with guide probe (3 x 240s exposure). Therefore, (B) and 
(C) are independent and taken under different observing 
modes, while (A) is a combination of (B) and (C). We 
note that the A520 data were taken with large dithering 
offsets of « 2f3, whereas the other cluster data were taken 
with a dithering offset of 1'. 

Our weak lensing analysis is based on co-added mosaics 
of multiple exposures (A), (B), and (C), whereas Mahdavi 
et al. (2007) performed shape measurements separately 
for each of seven Subaru exposures and combined individ- 
ual weak lensing catalogs to improve the accuracy for the 
shear estimates. In the present paper, the deepest data 
set (A) is used for our main analysis of A520, and the 
resulting mass map is compared with X-ray and optical- 
luminosity distributions of the cluster in Figure 11. Figure 
18, on the other hand, compares the mass reconstruction 
based on (B) with the X-ray and optical-luminosity dis- 
tributions. 



We show in Figure 19 the mass maps of A520 recon- 
structed from the above three imaging data sets. The left 
panel in Figure 19 shows the mas map from the deep- 
est data set (A), which is used for the main analysis of 
the present paper (see also Figure 11); the middle panel 
shows the mass map from (B) taken without guide probe; 
the right panel shows the mass map from (C) taken with 
guide probe. For all of the reconstruction, the Gaussian 
smoothing FWHM is taken to be lf25. The contours are 
spaced in units of Itr reconstruction error of each data set. 

The primary mass peak CI found from the data set (A) 
is also detected at a significance of 5.6cr from the data set 
(C) taken with guide probe. The angular position of CI 
agrees with that of the dark core discovered by Mahdavi 
et al. (2007) within Gaussian FWHM used for the map- 
making. However, this mass clump is not particularly 
pronounced in the mass map obtained from (B). 

The mass clump C3 of (A) is associated with a local con- 
centration of cluster member galaxies (see Figure 11). The 
mass clump C3 is found to be the first peak in the recon- 
structed mass map from (B). The mass structure around 
C3 is not strongly pronounced in the mass map based on 
(C), but the reconstruction shows a clear feature of the 
extending mass structure towards the southeast direction, 
which is seen in the optical luminosity map of the cluster 
(see the bottom-left panel of Figure 11). However, we note 
Mahdavi et al. (2007) found no significant mass struc- 
ture around C3 from their multi-telescope, multi-bandpass 
weak lensing analysis based on Subaru/Suprime-Cam Rc 
and i' and CFHT/MegaCam r' data. 

The mass clump C2 of (A) corresponds to the primary 
peak in the optical luminosity distribution of the cluster. 
The mass reconstruction based on (B) shows a local max- 
imum around C2, whereas no local maximum is found 
around C2 in the mass map from (C). 

The mass structures SW2, NEl, and NE2 are seen in 
all of the three reconstructions above a significance level 
of 4(T. The mass clump SWl is detected at a significance 
level of 3.5(7 from (A), and is marginally detected at « 3a 
significance from (B), but is not seen in the mass recon- 
struction from (C). 

Details of mass structures are different in a qualita- 
tive and quantitative way between the results from (B) 
and (C) taken without and with guide probe, respectively. 
This difference in detailed mass structures might be due 
to the different pointing modes and exposures between 
the data sets (B) and (C), which could affect the shape 
measurements of faint background galaxies. Furthermore, 
the large dithering offsets of 2. '3 could also make it diffi- 
cult to correct for the PSF anisotropy in stacked images, 
as discussed by Mahdavi et al. (2007). Careful follow-up 
imaging observations with guide probe and small dithering 
offsets (^ 1') will make it possible to clarify these discrep- 
ancies in detailed mass structures in A520. 
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Table 3. Color-Magnitude relation 



Cluster 


a 


1. 

b 


c 


T It 

-tJtot/ J^ohs 


Uj 


[^) 








A754 


-0.00999 


1.15662 


0.16 


1.00103 


A1750 


-0.02849 


1.65197 


0.12 


1.00268 


A1758 


-0.11661 


3.93112 


0.16 


1.04095 


A1914 


-0.04522 


2.13056 


0.10 


1.01238 


A2034 


-0.02424 


1.49218 


0.18 


1.00483 


A2142 


-0.02019 


1.51493 


0.11 


1.00302 


A520 


-0.0236 


1.33718 


0.12 


1.00531 



Note — Col. (1): Cluster name. Col. (2): 
Color-magnitude relation parameters of the form: 
jcolor — (a X mag + h)\ < c, where color = g' — Rc 
and mag = Rc for all targets, except color ~ V — i' 
and mag = i' for A520. Col. (3): Correction factor 
of galaxy luminosities due to the luminosity cutoff, 
itot/iobs - T{2-p)/r{2-p,Lnm/L*). 
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Fig. 1. The color-magnitude diagram of A2142: The early- 
type galaxies in clusters are well separated by their colors 
thanks to the redshifted 4000A break. 
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Table 1. Target Clusters and Subaru /'K-ra.y Observations 



Cluster 
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& Exposure Time 
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A754 


i?c = 360 s X 8 


5' = 180 s X 4 


0.75 


Chandra 


577 


A1750 


i?c = 360 s X 8 


5' = 180 s X 4 


0.67 


XMM 


0112240301 


A1758 


i?c = 360 s X 8 


g' = 180 s X 4 


0.69 


XMM 


0111160101 & 0142860201 


A1914 


i?c = 360 s X 8 


g' = 180 s X 4 


0.61 


Chandra 


3593 


A2034 


i?c = 360 s X 8 


g' = 180 s X 4 


0.63 


XMM 


0149880101 


A2142 


i?c = 240 s X 4 + 360 s X 3 


g' 120 s X 4 


0.55 


Chandra 


5005 


A520^ 


i' = 240 s X 7 


l/==450 s X 4+120 s X 3 


0.65 


Chandra 


4215 



Note — Col. (1): Cluster name. Col. (2): Band name of Subaru/Suprime-Cam and total exposure time in units of 
s. The former passband is used for the weak lensing analysis. Col. (3): Seeing FWHM in Rc or i' band in units of 
arcsec. Col. (4): X-ray satellite name of archival data. Col. (5): X-ray observation ID number. 

Data were collected at Subaru telescope and obtained from the SMOKA science archive. Four i'-band images were 
taken without guide probe on 19th October 2001, while the other three i' images were taken without guide probe on 
17th November 2001. 



Table 2. Cluster X-ray Properties 



Cluster 


z 


Type 


1 arcmin 


Components 


^avc 








(kpc /ifo^) 




(keV) 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


A754 


0.0542 


On-going 


63.1 




10.0 ±0.3^ 


A1750 


0.0860 


Binary 


96.7 


A1750C 


3.87±0.10 










A1750N 


2.84±0.12 


A1758 


0.2790 


Binary 


254.0 


A1758N 


8.2 ±0.4 










A1758S 


6.4t°i - 


A1914 


0.1712 


On-going 


174.9 




10.9±0.7 ^ 


A2034 


0.1130 


Cold Front 


123.2 




7.9 ±0.4 


A2142 


0.0909 


Cold Front 


101.7 




8.1 ±0.4 


A520 


0.1990 


On-going 


197.2 




7.1 ±0.9 



Note — Col. (1): Cluster name. Col. (2): Cluster redshift. Col. (3): Cluster merger type, classified by their X-ray 
features. 'Binary' represents the system consisting of two clusters separated by the order of Mpc scales. A 'Cold 
Front' cluster shows the ICM with the presence of a contact discontinuity. 'On-going' is a cluster whose X-ray surface 
brightness and temperature maps show irregular morphologies. Col. (4); Physical scale in kpc/Ziyo unit corresponding 
to 1' at the cluster redshift. Col. (5): Name of a component of binary cluster. Col. (6): Spatially averaged temperature 
from Govoni et al. (2004), Belsole et al. (2004), = David & Kempner (2004), ^ Kempner & Sarazin (2003) and ^ 
Markevitch et al. (2000). 
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Table 4. Weak Lensing Stellar Sample 



Cluster 


Uncorrected e* 






Residual 5e* 








N* 






el X 10^ 


el X 10^ 




X 10^ 


5 el X lO'' 


a{5e 


*) 




arcsec 


(1) 


(2) 






(3) 








(4) 


(5) 


A754 


-2.10 


-0.52 


1.77 X 10^^ 


+0.63 + 1.30 


+0.10 + 0.88 


4.25 X 10 


-3 


733 


0.440 


A1750 


+0.65 


-0.05 


1.73 X 10-2 


-0.46 + 1.54 


-0.11 + 0.92 


4.60 X 10" 


-3 


651 


0.384 


A1758 


+0.10 


+0.14 


1.77 X 10-2 


0.00 + 2.44 


+0.06 + 1.00 


4.20 X 10" 


-3 


255 


0.380 


A1914 


-0.85 


+0.50 


2.25 X 10-2 


+0.42 + 2.67 


-0.18 + 1.68 


6.11 X 10- 


-3 


375 


0.345 


A2034 


-1.03 


+2.13 


2.51 X 10-2 


+0.02 + 2.82 


-2.70 + 2.10 


7.52 X 10- 


-3 


458 


0.342 


A2142 


-1.23 


-2.77 


2.86 X 10-2 


+1.01 + 2.64 


+4.24+1.82 


8.39 X 10 


-3 


684 


0.294 


A520 


-2.40 


-2.43 


2.98 X 10-2 


+0.19 + 1.55 


+2.10 + 1.29 


7.92 X 10 


-3 


1537 


0.374 



Note — Col. (1): Cluster name. Col. (2): Mean and standard deviation of stellar ellipticities before PSF correction. 
Col. (3): Mean and standard deviation of stellar ellipticities after PSF correction. Col. (4): Number of stellar 
objects. Col. (5): Median stellar half-light radius. 



Table 5. Background Galaxy Sample and Weak Lensing Mass Reconstruction 



Cluster 


Magnitude range 


Tig 


^9 


FWHM 






ABmag 


arcmin 2 




arcmin 




(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


A754 


19 < i?c ;S 27 


37.3 


0.376 


1.67 


0.01739 


A1750 


19 < i?c ;S 27 


56.4 


0.405 


1.25 


0.02032 


A1758 


19 < i?c ^ 27 


46.4 


0.408 


1.25 


0.02253 


A1914 


19 < i?c ;S 27 


47.7 


0.411 


0.75 


0.03737 


A2034 


19 < i?c ;S 27 


52.4 


0.413 


1.17 


0.02301 


A2142 


19 < i?c ;S 27 


72.1 


0.395 


1.00 


0.02190 


A520 


20 < i' < 25.5 


37.6 


0.424 


1.25 


0.02605 



Note — Col. (1): Cluster name. Col. (2): Magnitude range of the background galaxy sample. Col. (3): Surface 
number density of background galaxies for map-making. Col. (4): RMS error for the shear estimate per galaxy. Col. 
(5): Gaussian FWHM in map-making. Col. (6): RMS noise in the reconstructed mass map. 



Table 6. Red Background Galaxy Sample 



Cluster 


Magnitude range 


Color 


Ug 




ABmag 




■ -2 

arcmm 


(1) 


(2) 


(3) 


(4) 


A754 


19 < i?c ;S 27 


1.14<.g'-i?c < 15.0 


6.61 


A1750 


19 < i?c ;S 27 


1.23<.g'-i?c < 15.0 


8.55 


A1758 


19 < i?c ^ 27 


1.20<.g'-i?c < 15.0 


7.90 


A1914 


19 < i?c ^ 27 


1.36<.g'-i?c < 15.0 


5.11 


A2034 


19 < i?c ;S 27 


1.20<.g'-i?c < 15.0 


7.31 


A2142 


19 < i?c ^ 27 


1.29<g'-i?c < 15.0 


7.55 


A520 


20 < i' < 25.5 


0.90 <y-i'< 15.0 


9.40 



Note — Col. (1): Cluster name. Col, (2): Magnitude range of the red background galaxy sample. Col. (3): Color 
range of the red background galaxy sample. Col. (4): Mean surface number density of the red background galaxies. 
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Table 7. Best-Fitting Mass Models 



Cluster 


SIS 






NFW 










ft) 


xVdof (dof) 






Cvir 


xVdof (dof) 






(km s^-'^) 




(arcmin) 


(/i-ilOi^Mo) 






(arcmin) 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


A754 


703.1 ±139.5 


0.32(3) 


(6.7,15.0) 


2.86 ±3.77 


4.97±5.26 


2.57(4) 


(1.5,5.4) 


A1750N 


633.9 ±80.2 


1.30(7) 


(1.0,7.0) 


2.90 ±2.08 


5.21 ±3.42 


1.11(6) 


(1.0,7.0) 


A1750C 


588.9 ±90.6 


0.91(7) 


(1.0,6.7) 


1.88 ±1.63 


5.52 ±4.89 


0.99(6) 


(1.0,6.7) 


A1758N 


672.5 ±113.3 


1.23(4) 


(1.7,5.0) 


5.26±5.70 


0.24 ±0.95 


1.01(5) 


(1.0,5.0) 


A1758S 


594.5 ±94.0 


3.51(7) 


(0.5,4.6) 


1.20±1.15 


3.91±6.51 


1.04(6) 


(0.5,4.6) 


A1914 


845.7 ±87.4 


1.08(7) 


(0.5,5.4) 


6.14±3.19 


4.13±2.79 


0.75(6) 


(1.0,7.0) 


A2034 


792.8 ±78.2 


1.46(10) 


(1.0,12.0) 


7.17±4.30 


3.74 ±2.25 


0.84(9) 


(1.0,12.0) 


A2142 


940.4 ±62.6 


2.03(10) 


(1.0,15.0) 


11.95±5.24 


4.32 ±1.70 


0.83(9) 


(1.0,15.0) 


A520 


960.4 ±94.0 


0.70(4) 


(5.1,15.0) 


6.14±2.38 


2.90 ±1.82 


0.68(7) 


(1.7,15.0) 



Note — Col. (1): Cluster name. Col. (2): Best-fitting parameter value for the SIS ID velocity dispersion. Col. 
(3): Reduced for the best-fitting SIS model, and the dcgrccs-of- freedom in parenthesis. Col. (4): Radial range 
of the Cc-statistic measurements used for the SIS model fitting. Col. (5): Best-fitting parameter value for the NFW 
virial mass. Col. (6): Best-fitting parameter value for the NFW concentration parameter. Col. (7): Reduced for 
the best- fitting NFW model, and the degrees-of-freedom in parenthesis. Col. (8): Radial range of the Cj,-statistic 
measurements used for the NFW model fitting. 



Table 8. Global Cluster Properties 



Cluster 


(/i-ilOi^M©) 


L{<9m) 


MjLi< 9m) 




(^inn, ^out) 




(HMq/Lq) 


(arcmin) 


(arcmin) 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


A754 


2.73±1.70 


0.291^:^^ 


932+^^5 


15.0 


(15.0,20.0) 


A1750N 


1.49 ±0.65 




660+2^° 


6.96 


(16.0,21.0) 


A1750C 


1.26±0.60 


"■r' -0.01 


7f;n-|-366 


6.65 


(15.0,20.0) 


A1758N 


4.89 ±1.54 


9 -1 7-1-0.13 
^■^ ' -0.02 




5.08 


(15.0,20.0) 


A1758S 


6.06±1.10 


i.2ot°:J« 


50611^6 


4.63 


(12.0,17.0) 


A1914 


4.10±1.55 


1.081°:- 




6.96 


(15.0,20.0) 


A2034 


3.90±1.81 


-■^■0^-0.05 


qoo-f-lTS 


12.0 


(12.0,17.0) 


A2142 


3.97 ±2.04 


1 09+1-37 

-'^•o-'-o.is 


365_274 


15.0 


(15.0,20.0) 


A520 


5.18±4.19 


o An+1.35 


1 r;9+123 
-•^^^-130 


15.0 


(15.0,20.0) 



Note — Col. (1): Cluster name. Col. (2): Cluster mass A/^(< 9m) enclosed within 9m from aperture densitometry. 
Col. (3): Cluster luminosity enclosed within 9m in i?c-band for all targets, except in i'-band for A520. Col. 
(5): Cluster mass-to-light ratio within 9m- Col. (6): Aperture radius in units of arcmin. Col. (8): Inner and 
outer boundaries of the background annulus for the cluster mass and luminosity measurements, given in units of arcmin. 
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Fig. 2. The lower panel shows the inverse of the critical surface mass density of gravitational lensing, ^^j.^ {z^,Zs), as a function of 
lens redshift for three different source redshifts, Zs = 0.8,1.0,1.2 {dashed, solid, and dotted- dashed, respectively), demonstrating 
the geometric scaling of gravitational lensing signal. The top panel shows the relative lensing strength Tj^^^ (z^, Zs) /l^cr^ (z^i, Zs = 1.0) 
as a function of lens redshift 2jj normalized with respect to the source at z^ = 1. For lensing clusters at low redshifts z^, S^r depends 
very weakly on the background redshift Zs , so that the uncertainty in Zs of background galaxies is less important in the lensing-based 
cluster mass determination. 
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Fig. 3. Stellar ellipticity distributions before and after the PSF anisotropy correction for individual cluster targets. For each 
cluster the left panel shows the raw ellipticity components (e^jCj) of stellar objects, and the right panel shows the residual ellipticity 
components (5e*,(5e2) after the PSF anisotropy correction. For A520 two different imaging data sets are shown: A520 (A) is based 
on a co-added mosaic of seven i' images (7 X 240s) taken both with and without guide probe. A520 (B) is based on a co-added 
mosaic of four i' images (4 x 240s) taken without guide probe. 
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Fig. 4. The distortion field of stellar ellipticities before and after the PSF anisotropy correction for individual cluster targets. For 
each cluster the left panel shows the raw ellipticity field of stellar objects, and the right panel shows the residual ellipticity field 
after the PSF anisotropy correction. The orientation of the sticks indicates the position angle of the major axis of stellar ellipticity, 
whereas the length is proportional to the modulus of stellar ellipticity. A stick with the length of 10% ellipticity is indicated in the 
top left of the left panel for each target. For A520 two different imaging data sets are shown: A520 (A) is based on a co-added 
mosaic of seven i' images (7 X 240s) taken both with and without guide probe. A520 (B) is based on a co-added mosaic of four i' 
images (4 X 240s) taken without guide probe. 
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Fig. 5. A754 — Top-left: Subaru _Rc-band image of the central ~ 25' X 13' cluster region. Overlayed are contours of the lensing 
K- field reconstructed from weak shear data. The contours are spaced in units of Icr reconstruction error (see Table 5). The Gaussian 
FWHM used for the mass reconstruction is 1.'67. Top-right: adaptively-smoothed Chandra X-ray contours (0.7 — 7.0 keV) overlaid on 
the same iJc-band image. Bottom-left: Cluster luminosity density distribution in iJc-band smoothed to the same angular resolution 
of the mass map. Overlayed are the same mass contours as in the top-left panel. Bottom-right: The same mass contours overlayed 
on the adaptively-smoothed Chandra X-ray image (0.7 — 7.0 keV). 
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Fig. 6. A1750 — Top-left: Subaru iJc-band image of the central ~ 21', xl8' cluster region. Overlayed are contours of the lensing 
K-field reconstructed from weak shear data. The contours are spaced in units of Icr reconstruction error (see Table 5). The Gaussian 
FWHM used for the mass reconstruction is 1^25. Top-right: adaptively-smoothed XMM-Newton X-ray contours (0.5 — 7.0 keV) 
overlaid on the same Rc-hand image. Bottom-left: Cluster luminosity density distribution in i?c-band smoothed to the same angular 
resolution of the mass map. Overlayed are the same mass contours as in the top-left panel. Bottom-right: The same mass contours 
overlayed on the adaptively-smoothed XMM-Newton X-ray image (0.5 — 7.0 keV). 
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Fig. 7. A1758 — Top- left: Subaru i?c-band image of the central ~ 26', xl9' cluster region. Overlayed are contours of the lensing 
K- field reconstructed from weak shear data. The contours are spaced in units of la reconstruction error (see Table 5). The Gaussian 
FWHM used for the mass reconstruction is lf25. Top-right: adaptively-smoothed XMM-Ncwton X-ray contours (0.5 — 7.0 kcV) 
overlaid on the same i?c-band image. Bottom-left: Cluster luminosity density distribution in iJc-band smoothed to the same angular 
resolution of the mass map. Overlayed arc the same mass contours as in the top-left panel. Bottom-right: The same mass contours 
overlayed on the adaptively-smoothed XMM-Newton X-ray image (0.5 — 7.0 keV). 
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Fig. 8. A1914 — Top-left: Subaru i?c-band image of the central ~ 16', xl5.'5 cluster region. Overlayed are contours of the lensing 
K- field reconstructed from weak shear data. The contours are spaced in units of Icr reconstruction error (see Table 5). The Gaussian 
FWHM used for the mass reconstruction is Of 75. Top-right: adaptively-smoothed Chandra X-ray contours (0.7 — 7.0 kcV) overlaid on 
the same iJc-band image. Bottom-loft: Cluster luminosity density distribution in i?,c-band smoothed to the same angular resolution 
of the mass map. Overlayed arc the same mass contours as in the top-left panel. Bottom-right: The same mass contours overlayed 
on the adaptively-smoothed Chandra X-ray image (0.7 — 7.0 keV). 
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Fig. 9. A2034 — Top-left: Subaru Kc-band image of the central ~ 18f5, xl8'.5 cluster region. Overlayed are contours of the lensing 
K- field reconstructed from weak shear data. The contours are spaced in units of la reconstruction error (see Table 5). The Gaussian 
FWHM used for the mass reconstruction is 1'.17. Top-right: adaptively-smoothed XMM-Newton X-ray contours (0.5 — 7.0 keV) 
overlaid on the same i?c-band image. Bottom-left: Cluster luminosity density distribution in iJc-band smoothed to the same angular 
resolution of the mass map. Overlayed are the same mass contours as in the top-left panel. Bottom-right: The same mass contours 
overlayed on the adaptively-smoothed XMM-Newton X-ray image (0.5 — 7.0 keV). 
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Fig. 10. A2142 — Top-loft: Subaru /?c-band image of the central ~ 25', x23f5 cluster region. Overlayed are contours of the lensing 
re- field reconstructed from weak shear data. The contours are spaced in units of Itr reconstruction error (see Table 5). The Gaussian 
FWHM used for the mass reconstruction is IfOO. Top-right: adaptively-smoothed Chandra X-ray contours (0.7 — 7.0 keV) overlaid on 
the same iic-band image. Bottom-left: Cluster luminosity density distribution in iJc-band smoothed to the same angular resolution 
of the mass map. Overlayed are the same mass contours as in the top-left panel. Bottom-right: The same mass contours overlayed 
on the adaptively-smoothed Chandra X-ray image (0.7 — 7.0 keV). 
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Fig. 11. A520 — Top-left: Subaru i'-band image of the central ~ 12', xl2' cluster region. Overlayed are contours of the lensing 
K- field reconstructed from weak shear data. The contours are spaced in units of Ifr reconstruction error (see Table 5). The Gaussian 
FWHM used for the mass reconstruction is 1.'25. Top-right: adaptively-smoothed Chandra X-ray contours (0.7 — 7.0 keV) overlaid 
on the same j'-band image. Bottom-left: Cluster luminosity density distribution in i'-band smoothed to the same angular resolution 
of the mass map. Overlayed are the same mass contours as in the top-left panel. Bottom-right: The same mass contours overlayed 
on the adaptively-smoothed Chandra X-ray image (0.7 — 7.0 keV). 
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Fig. 12. Radial profiles of the projected mass Mq{6) for a sample of seven merging clusters of galaxies as measured by the weak 
lensing ^-statistic. The error bars are correlated. Also plotted are the best-fitting NFW and SIS profiles. 
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Fig. 13. Comparison between the best-fitting SIS temperature Tsis from the (^c-statistic measurement vs. the averaged X-ray 
temperature Tx for a sample of nearby merging clusters. The dashed line indicates Tsjs = Tx- For all clusters except A1750N, the 
observed X-ray temperature is significantly higher than the expected virial temperature of the cluster mass. 
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Fig. 14. The temperature ratio Tx/Tgig v.s. the best-fitting SIS temperature Tgjg. The dashed line indicates Tgjg = Tx- 
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Fig. 15. Subaru Rc-hand image of the central 8.' X 8.' cluster region of A1758N. Overlayed are contours of the reconstructed projected 
mass distribution of the cluster. The insets zoom in on tangential arc candidates identified based on visual inspection. The panels 
A (0.'4 X 0.'4) and B {0'.2 X 0'.2) are zoom in of blue tangential arcs around the SE mass clump. The panel C (0.'2 X 0.'2) is a zoom 
in view of blue tangential arcs around the mass clump C. The panels D (0.'3 X 0.'3) and E (0.'3 X 0'.3) arc zoom in of tangential arcs 
associated with cluster galaxy concentrations. 




Fig. 16. Subaru _Rc-band image of the central 5.' X 5.' cluster region of A1914. Overlayed arc contours of the reconstructed projected 
mass distribution of the cluster. The inset panel (Of4 X 0.'4) shows a zoom in view of two blue tangential arcs around BCG2. 
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Fig. 17. Subaru i?c-band image of the central 8^ X 8'. cluster region of A2142. Ovcrlaycd arc contours of the reconstructed projected 
mass distribution of the cluster. The insets zoom in on tangential arc candidates identified based on visual inspection. The panel A 
shows a zoom in view (0f4 X 0.'4) of a tangential arc. The panel B shows a zoom in view (0.'4 X 0.'4) of a tangential arc. The panel 
C sows a zoom in view (0.'8 X 0'.8) around BCG2, where no arc-like image is found by visual inspection. 
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Fig. 18. A520 based on 4 X 240s j'-band images taken without guide star — Top-left: Subaru i'-band image of the central ~ 11', xll' 
cluster region. Overlayed are contours of the lensing K-ficld reconstructed from weak shear data. The contours are spaced in units of 
1(T reconstruction error (see Table 5). The Gaussian FWHM used for the mass reconstruction is 1^25. Top-right: adaptively-smoothed 
Chandra X-ray contours (0.7 — 7.0 keV) overlaid on the same i'-band image. Bottom-left: Cluster luminosity density distribution in 
i'-band smoothed to the same angular resolution of the mass map. Overlayed are the same mass contours as in the top-left panel. 
Bottom-right: The same mass contours overlayed on the adaptively-smoothed Chandra X-ray image (0.7 — 7.0 keV). 
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Fig. 19. Comparisons of A520 weak Icnsing analyses of different data sets. Left (A) : The lensing K-field of (12.' X 12.'), reconstructed 
from seven i' images with and without AG (acquisition and guide) probe (7 X 240s exposure), which was also used by Mahdavi et 
al. (2007). The Gaussian FWHM used for the mass reconstruction is 1.'25. The contours are spaced in units of la reconstruction 
error. Middle (B) : The same figure as left panel, reconstructed from four i' images taken without guide probe (4 X 240s exposure). 
Right (C) : The same figure as left panel, reconstructed from three i' images taken with guide probe (3 X 240s exposure). 
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